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AN ANALYTICAL INVESTIGATION OF WING- JET INTERACTION 


by C. Edward Lan 
The University of Kansas 
SUMMARY 

The aerodynamic interaction between the wing and an invlscid jet 
with Mach number nonuniformity is formulated by using a two-vortex- 
sheet model for the jet. One of the vortex sheets accounts for the 
Induced jet flow and the other the induced outer flow. No additional 
source distribution is needed for the jet at an angle of attack. The 
above problem is solved by satisfying the jet and wing tangency and the 
jet pressure-continuity conditions and using a "quasi vortex lattice 
method" for computing the induced flow field. The latter method is 
derived through theoretical consideration by properly accounting for 
singularities present in the equations and possesses the same simplicity 
and generality as the conventional vortex lattice method but has a better 
rate of numerical convergence. The resulting system of algebraic equa- 
tions is solved by Purcell’s vector method. The numerical formulation 
is first applied to the wing-slipstream interaction problem. Results 
for one centered-jet configuration are compared with those predicted by 
some existing theories. The present method is shown to have the same 
accuracy as Rethorst’s. Good results have also been obtained in lift, 
induced drag and pitching moment coefficients when compared with available 
experiments. If the slipstream Mach number is taken to be greater than 
the free stream value, the loading is shown to be increased, but the incre- 
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merit Is much less than expected from Prandtl-Glauert transformation 
applied locally. Explanation is given for this characteristic. The 
method is finally applied to the upper-surface-blowing STOL configu- 
ration. The agreement of the present prediction with measurement by 
Phelps, et al in lift Increment is good. Various characteristics of 
the aerodynamic interaction have been pointed out and discussed. 



1. Introduction 


In the design of early subsonic turbojet aircraft, the inter- 
ference effects of airframe and engines were not regarded as major 
problems. However, with the advent of high bypass ratio turbofan 
engines, significant interference effects have been recognized. In 
fact, the integration of the propulsion system with the airframe has 
been regarded as one of the major causes which create the discrepancy 
in the predicted performance with the delivered one. Schumacher and 
Trent [1] analyzed the performance of some existing fighters and found 
that the increased engine size is definitely to reduce the cruise ef- 
ficiency. Attempt has been made recently through wind-tunnel testing 
with propulsion simulator to identify the complex mechanisms of air- 
frame/engine interference [2,3,4]. Measurements of two-dimensional 
pressure distributions due to jet effects for the conventional under- 
wing engine installation were made at high subsonic Mach numbers [5]. 

It was found that the jet affects only the lower surface pressure 
distribution and that the jet effect can not be adequately simulated 
by an equivalent solid body. Bagley further suggested that in order 
to satisfy the ground clearance requirement for large turbofan engines, 
the engines may be installed above the wing [6]. He found in the 
experiment on two-dimensional wings \/ith the jet pipe flattened to a 
rectangular shape near the exit that both lift and lift-dependent drag 
are increased for the above-wing engine installation. Similar engine 



configurations for the purpose of STOL operations have also been tested 
at NASA Langley Research Center [7]. In these latter engine installa- 
tions, it is obvious that accounting of the jet interaction becomes a 

necessity in predicting the wing aerodynamic characteristics. 

of low 

In general, the jet effects^may be divided into (1) the entrain- 
ment effect with external fluid being sucked into the jet, (2) the 
displacement effect, the external streamlines being displaced due to 
the presence of the jet, and (3) the interaction effect, in which the 
jet flow that is affected by the wing flow field influences in turn the 
wing flow. No general analytical methods of incorporating these effects 
in the aerodynamic computation are available. Early analytical or semi- 
empirical techniques to assess the gross jet effects on the performance 
of aircraft are discussed in Refs, [8] and [9]. More recently, Men- 
denhall and his co-workers investigated the STOL wing characteristics 
with externally-blown jet-augmented flaps [10]. They postulated that 
the jet effects can be accounted for by a vortex ring distribution on 
the jet boundary. The jet boundary assumed by them is in fact the 
boundary at which the longitudinal velocities of the jet and the external 
stream are equal, as has been defined and used in the turbulent jet 
theory. The variation of the vortex ring strengths downstream was 
obtained by assuming that the product of the strength and the ring cir- 
cumference remains constant. The result was that the conservation of 
jet momentum was violated. In addition, no interference effect, as 
defined above, has been included. A more rigorous approach in eval- 
uating the interaction effect was employed by Shollenberger [11 a & b]. 
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The location and strength of flow singularities were determined through 
iteration. However, the jet was assumed to be inviscid and incompres- 
sible. It is not known how important the jet distortion effects are 
on the wing loading. 

Historically, interest in the interaction of two flows with dif- 
fenent total energy levels has been centered on the problem of pro- 
peller-wing interaction. To gain some insight on the analytical 
techniques used in this problem, it is appropriate to review some of 
the developments. Koning [12] explained in great detail in 1935 how 
the slipstream boundary conditions of flow tangency and pressure con- 
tinuity can be satisfied by introducing additional flow singularities. 
Ferrari extended the formulation to the compressible flow through 
Prandtl-Glauert equation and in the spirit of the classical lifting 
line theory [13]. Since then, efforts have been directed toward the 
refinement of the theory by using the Weissinger-type lifting surface 
method. A main feature of these efforts is the employment of the 
image technique to find the additional flow singularities necessary 
to satisfy the slipstream boundary conditions. A representative of 
this type of work is Rethorst’s paper [14] which is concerned with a 
circular slipstream. Jameson [15] extended Rethorst*s method to treat 
slipstreams of rectangular cross sections, Ribner and Ellis [16] 
departed from the tradition of the image technique by presenting a 
formulation applicable to undeflected multi-slipstream configurations 
by using surface doublet distributions on the slipstreams. Even though 
the formulation is general enough for lifting-surface type computation. 
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computer implementation of the general case was not successful. Later, 
a "multiple lifting line approximation" was employed [17], in which both 
the wing and the slipstream were divided into streamwise segments. In 
each segment, a lifting line or lifting ring was placed at the quarter 
chord. However, the boundary conditions were satisfied only at the 
three-quarter chord of the rearmost segment. The relative values of 
the strengths of lifting lines or lifting rings were determined by 
matching the upwash. No accurate pressure distribution is obtainable 
by their program. Levinsky and others [18] extended the theory of 
Ribner and Ellis [16] to the case with deflected slipstreams. They 
represented the slipstreams by ring vortices and multipole distribu- 
tion on the slipstream axis. The slipstream boundary conditions are 
satisfied in the far field in the linearized approximation through the 
use of surface distribution of discrete horseshoe vortices. The wing 
aerodynamic characteristics were computed by the discrete-type Weissinger’s 
method. All wing and slipstream trailing vortices were assumed straight 
and inclined downward in the extended slipstream direction. No inter- 
ference between slipstreams was assumed. Again, the wing pressure 
distribution was not predicted. More recently, Kleins tein and Liu [19] 
showed that the incorporation of more accurate airfoil data in the 
lifting line theory for the wing-propeller interaction is important. 

A systematic procedure through the method of matched asymptotic expan- 
sions for the interaction of wing and multipropellers was later presented 
for high aspect-ratio wings [20]. 
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From the above review of literatures in the propeller wing inter- 
action, it is clear that the available methods, except Ref. [11b], are 
not adequate for predicting the wing pressure distribution and the 
spanwise induced drag distribution for arbitrary planforms without 
restrictions on aspect ratios and sweep angles. In addition, no method 
to account for the differences in the jet and external Mach numbers 
seems available. In the present investigation, the following improve- 
ments over the existing techniques will be made; (1) A quasi vortex 
lattice model for the wing will be used. This is a quasi continuous 
vortex model and represents an improvement over the conventional vor- 
tex lattice method. The wing pressure distribution and the spanwise 
induced drag distribution can be computed for arbitrary wing planforms. 
(2) The differences in the jet and external Mach numbers are accounted 
for through the use of two vortex sheets on the jet boundary. (3) 

The jet surface boundary conditions are satisfied everywhere on the 
jet surface, not just in the wing region as done in Ref. [17], The 
restrictions are; (1) The boundary conditions are satisfied only in 
the linearized approximation. Thus, the jet boundary conditions 
are satisfied on the jet surface and the tangency and pressure condi- 
tions are linearized. (2) The wing thickness is not included. (3) 

The Prandtl-Glauert equation is used. 
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2, List of Symbols 


a 



A 


b 



(xi-x)l + (yi-y)j + (zi-z)k, m(ft) 

(xi-x)i + 3(yi-y)j + 3(zi-z)k, m(ft) 
aspect ratio 

jet cross-sectional area, m^(ft^) 
span, m(ft) 

(x 2 ~x)i + (y 2 -y)j + (z 2 -z)k, m(ft) 

(x 2 -x)i + 3 (y 2 “y)j + 6 (z 2 ~z)k, m(ft) 
local chord length, m(ft) 
reference chord length, m(ft) 

leading edge suction parameter. See Eqs. (4.25), (4.26) 

sectional induced drag coefficient 
total induced drag coefficient 
drag coefficient due to Coanda effect 
sectional lift coefficient 

sectional lift coefficient based on average dynamic pressure; 

is equal to 2 Cj^p ^/(1 + y^) 

total lift coefficient 

lift coefficient due to Coanda effect 

sectional pitching moment coefficient about the y-axis 

total pitching moment coefficient about the y-axis 

pressure coefficient 

sectional leading-edge thrust coefficient 
total leading-edge thrust; or propeller thrust coefficient 
equal to thrust divided by the disc area and the free stream 
dynamic pressure 
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Jet thrust coefficient, T /q S 

r 

diameter of slipstream 
drag due to Coanda effect, N(lb) 
unit tangential vector 
a vector defined by Eq. (4.34) 

distance of the jet lower surface to the wing surface, m(ft) 

propeller advance ratio 

(x2"Xi)i + (y2“yi)j + (Z2^zi)k, m(ft) 

(^ 2 “Xi)i + B(y 2 -yi)j + e(z 2 -zi)k, m(ft) 

total lift, N(lb) 

lift due to Coanda effect, N(lb) 

Mach number 

a vector defined by Eq. (4.33) 
unit normal vector 

normal velocity influence coefficient matrix 
number of chordwise integration points 
static pressure, N/m^ (Ib/ft^) 
induced velocity vector, m/sec (ft/sec) 
dynamic pressure, N/m^ (Ib/ft^) 

A A 

xi + yj + zk, m(ft) 

(x-x')^ + B^(y-y')^ + 3^(z-z')^> m^(ft) 
radius of slipstream, m(ft) 
streamwise coordinate, m(ft) 

tangential velocity influence coefficient matrix 
wing area, m^(ft^) 
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z^(x) 

a 

a = 


propeller thrust coefficient based on q . See Eq. (5.12) 

s 

jet thrust, N(lb) 
jet thickness, m(ft) 

free stream velocity vector, m/sec (ft/sec) 
nondlmensional perturbed velocity in the x- direction 
unperturbed velocity vector, m/sec (ft/sec) 
perturbed velocity vector, m/sec (ft/sec) 
slipstream rotational velocity, m/sec (ft/sec) 
nondimens ional perturbed velocity in the 2 - direction 
wing-fixed rectangular coordinates with positive x- axis 
along the axis of symmetry pointing downstream, positive y- 
axis pointing to the right and positive z- axis pointing 
upward, m(ft) 
camber, m(ft) 
angle of attack, deg, 

(1-M^)^ 


612,621 diffraction coefficient. See Appendix B 
6 ^ flap deflection, deg 

y nondlmensional vortex density 

r discrete vortex strength 

K nondlmensional streamwise coordinate. See Eq. (4.36) 

^12»^21 reflection coefficient. See Appendix B 

A sweep angle, deg 

W = or V^/V^ 

W- = (5.S^)/(?3.e^) 
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dimensional perturbation velocity potential, m^/sec (ft^/sec) 
5 nondimensional perturbation velocity potential. See Eqs. 

(4.1), (4.2) 

p density, kg/m^ (slugs/ft^) 

nondimensional additional perturbation velocity potential. 

See Eqs. (4. 8) , (4. 9) 

0 Brenckmann’ s slipstream strength parameter equal to 2T 

c 

e angular coordinate. See Eq. (4.36) 

Subscripts 

1 first endpoint of a vortex element 

2. second endpoint of a vortex element 

a additional 

cp center of pressure 

j jet flow 

leading edge 
t trailing edge 

o outer flow 

s slipstream or jet region 

st strip theory 

w wing 
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3. Boundary Conditions 

The base flow, or the flow field before the wing is introduced, 
includes in general the uniform free stream, the jet flow and the jet 
entrained flow due to turbulent mixing. In order to simplify Che 
analysis, the jet entrained flow will be ignored in the present analysis. 
As the wing is introduced into the base flow, perturbations of the flow 
field occur. The calculation of these perturbations represent one 
of the main objectives in the present investigation. 

The perturbed flow field must be calculated by satisfying the 
boundary conditions. The boundary conditions on the jet surface 
require that the surface be a stream surface and the pressure there 
be continuous. The stream surface condition is satisfied if the 
slopes of streamlines on both sides of the jet surface are the same, 
i.e. 

K-{u + vi) + 

where e^ is a unit vector tangent to the jet path axis and n is a 
unit normal vector at the boundary of the jet cross section normal 

^ r, V 

to the jet path. See Figure 3.1. U is the freestream vector. V 

s 

is the jet flow velocity vector and xj and iT. are the perturbed velo- 
cities outside and inside the jet, respectively. 
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A-A 

Figure (3»1) Definition of Jet Coordinate System 


The pressure continuity condition can be found by using the 
Bernoulli’s equation: 


+ 2 ) ^tI Y- j 2 ‘ 


Using the isentropic condition: 

5'^rfef 

Eq. (3.2) can be reduced to 

= i + ^4^,[ivvr- i^tP3 

For the outer flow, let Pj-^£ = Pq , the pressure of the outer flow at 
the jet exit and Vj-^£ =» U . Then = U + V , For the jet flow, let 
Pref ^ Pj > 3^^ pressure at the jet exit, and ^ ^ 

Since p^ = pj , the continuity of the static pressure on the jet surface 


Implies that 


l'+-¥ff,citr-iv. »5, 

0 ’‘j 




If it is further assumed that , Eq. (3,5) can be reduced to 

C IvJ"- I V. f = ?.[ Ml"- \ u+ vj"3 

or, 

I vi - ( + Vj)l • C K + ( ■»■ Vj )) = 

§„ [ u - ( u + v„)j .[ n + ( u + u, )j 





Since the velocity difference across the surface can be represented 


by a vorticity distribution having the strength equal to the magnitude 
of the velocity difference and a direction normal to the velocity dif- 
ference, Eq. (3,6) can be interpreted as the continuity of Joukowsky 
force across the jet surface* Eq, (3,6) can also be recast to the 
following form; 

C Vs + = X [ U + -^ 3 

To simplify Eqs . (3,1) and (3,7) , let 

Uj. = ( U. eJCs -+ 

\ = v4>„ 

Substitution of these relations into Eq. (3.1) gives 

Cu-T> + v+.-Alf Vs-e. + = 

Cu- t Vs « + (o-e^K 

For the purpose of linearizing the last equation, it is assumed that 

« Vj -es C^-io) 

U • 6s (3.11) 

It follows that 

( Vs‘ - ( ij . 4)t ^ 

Then 


^3-k 

(3.7) 

gives 

■n)+ Ti- 74 - 5 } 

(3.6) 

(3.3) 


U-Ti ( 1-^0 -yw Vi-n ■f' ^ VPr 

where(a.^^'3(^. U-n has been used. Eq. (3,14) 

same as Eq, (8) of Ref. (18), Part I, if Vg- n = 0. 


C^. 12) 

(»,I3) 
(i. I4-) 


is exactly the 


or, 


Similarly, Eq, (3.7) can be written 

Sj-t V, + i(u. eje, + f (u-;,)e,+ V<j>^V= 

t [ O + i7'j>,3- vf. 

Sj [ ( u ej)( V,-e^) + 7 ( u-ej)’’ + ( U - e, t = 

S. ( U. 





etc. , it is 


Noting that ^ e.^) t (V^' V 


found that 


S,. ( V,- e,)^ + C V,-v.)% = S. t SAu.-^f^- 


■» n 


S-[(U-4H v^.ej +i-(‘J ezf+(u-eJv4>^-e2] C^'S) 


For small jet deflections, U.e^, U*n , V^-e^and \^‘n may be assumed 
to be of the same order as the perturbations. In this special case, 
Eq.(3.15) can be simplified to 


SAv,-4^)^ s 

1 S 




(3lfe) 


The wing boundary conditions will now be considered. Assume that 
the wing surface can be described by the following relation: 

^ ^ . y , 2;c;) - 2^ - C X , y ) — o I T) 

The unit normal vector is then 


yh ^ ^ i ^ 

Kill |yh| 


(3.18) 


The tangency condition requires that 



If the wing is not immersed in the jet, Eq,(3.19) can be written 





Since both tJ.| and are small, they may be neglected. 

Keeping only the first-order terms, it is found that 




(outside the jet) 


( 3 . 20 ) 


If the wing is inside the jet, the boundary condition Eq.(3.19) 
becomes 


C Vs + V. ) • ( 




From Eq.(3.8), it can be shown that 


(v^ + (,u-e^)e^ + - Ty 5' + = 0 


Or, approximately, 


DZ ^ DK 


- u- e^(.e, )V) - V^- 


^ -U-e,(v^)-Vs-fe 

(inside the jet) 
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4. Method of Flow Singularities 

Due to the presence of the wing, disturbances are created not 
only in the external flow, but also in the jet flow, which in turn 
affect the wing flow. Employing the technical terms associated with 
wave propagation, the interaction may be explained as the reflection 
and diffraction at the jet surface of the wing disturbances. The 
reflected disturbances in turn influence the wing flow and the dif- 
fracted ones perturb the jet flow. This principle is explained in 
more detail and applied to a two-dimensional idealized case in Appendix 
B. It is found in Appendix B that reflected disturbances account for 
most of the effects of wing- jet interaction for moderate or large 
jet thickness. In general, it may be stated that reflection of dis- 
turbances from the jet surface into the jet region will decrease the 
loading of that region from the uniform-flow value at the corresponding 
jet velocity. The contrary is true for the outer flow region. For 
the three-dimensional case with different Mach numbers inside and out- 
side the jet, quantitative prediction of these additional disturbances 
may be found by introducing two sheets of jet surface vortices, one 
being applicable in the outer region with Mq and one applicable in the 
jet region with M j , in such a way that boundary conditions Eqs.(3.14) 
and (3.16) are satisfied. The latter can be done in accordance with 
Koning*s method of introducing additional singularities. 

It is convenient to introduce nondimensional perturbation poten- 
tials 4’* and 4*^ such that 



(4.1) 


Eqs.(3.14), (3.16), (3.20) and (3.21) become 


^ H- U- n ( I 

-^y\ ~ (j-e^ -^T) 5. 




■Ti s 


- T(,--)"^ + TX {u %) 




;;> 

•^s 


(1^) 




U-e, 


(on wing surface) 


24l 



2^- 


( U-eJCe^-^) 




V^* e* 

(on wing surface) 




(^.O 


where 

T= S./S. 

In both Eqs.(4.3) and (4.4), it has been assumed that Vg’Cg is constant 
and Vg’n = 0. If it is further assumed that "2> ( Ci ♦ )/i n and ~b {, O ' j-cs> 

are small so that the corresponding terms are of the second order, 
then Eqs.(4.3) and (4.4) can be simplified to be 



2^ 


2t 

If n 




u ' n ( I -->*0 
U • Sj 






== T(A)’^ 


Ih 


(•t.4 ck) 


Let and ^l^.yjCMjv)be the nondimensional potential functions for 

the wing alone cases in the uniform flow xi^th the corresponding Mach 
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numbers of the outer flow and the jet flow, respectively. Let 



(4.S) 


C4f) 


It follows that and are the additional perturbation potentials 
for regions outside and inside the jet, respectively. Using Eqs.(4.8) 
and (4.9), the above boundary conditions become 


On jet J 

j ■» >1 

>n T* n 

CMP + -/;-^cM.) + ^ 

i 

A 

Ss 

(t.lo) 

surface 

'T* S 

'^>S 


(MO 

(4'l) 




•J A 1 

u ■ ft 



On wing ^ 


D ' 




surface 

'•>4i 



'V z. 


C U-ej.)Cea' /,) V.^ • ^ _ 




V, e, . 




In applications, Eqs . (4 ,10)-(4 .13) can be reduced to different forms 
for different problems. 


At this point, a two-vortex-sheet model for the jet surface will 
be introduced. The reasons for introducing two vortex sheets, instead 
of one, are discussed below. 

(1) If Ribner^s formulation [16] was followed, a single vortex 
sheet and a source sheet would have been introduced to account for 
the jumps in tangential velocities Eq. (4.11) and in normal velocicies 
Eq,(4.10), respectively. However, theoretically, these flow singulari- 
ties should be associated with the mean flow field which depends on 
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the unknown flow singularities. Thus, additional linearization as*“ 
sumption must be made to the effect that the mean flow field can be 
approximated by the unperturbed one. Furthermore, an additional 
different type of integral for the source distribution must be handled. 
The main advantage of Ribner^s formulation is that the number of equa- 
tions can be reduced by one. On the other hand, the present formu- 
lation will solve the problem defined in Eqs. (4,10)-(4.13) without 
additional linearization assumption. Besides, no source integral has 
to be introduced, because the jet and the outer perturbed flow fields 
will be solved separately by using the vortex distribution and the 
solutions matched according to Eqs, (4,10) and (4.11) at the jet surface. 
This is the same technique used by Rethorst [14], The disadvantage 
is that the number of unknowns to be determined would be increased. 

This disadvantage can be minimized by using an efficient technique for 
solving the simultaneous equations (4. 10)-(4.13) , to be discussed 
later in this section. 

(2) If then no existing methods are applicable. It is 

because if a flow singularity is introduced at a certain position on 
the jet surface, it will induce different normal and tangential velo- 
cities at all other points on the jet surface as they are approached 
from inside and outside the jet region. Therefore, the source and vor- 
tex distributions are interdependent. The number of unknowns can not 
be reduced. In this case, the present formulation is advantageous, 
because it is not necessary to introduce the source distribution. 


The present method can also be extended to the case where Mj > 1 and 

Mo < 1. 

Now define [Sjy] as the s-influence matrix for the jet due to 
the wing, i.e., the matrix elements being the s-perturbed velocities 
on the jet due to unit vortex distribution on the wing. Similar def- 
initions are applicable to the matrices [Nj^] , , [N^j] > etc. 

where N denotes the normal components. Using "y’* with appropriate 
subscripts to denote the unknown vortex strength, additional perturbed 
velocities on the jet surface can be written as 



(4:l<t) 


(4. fs) 


(4:IG) 




where the subscripts and superscripts ”o" and "j" denote, respectively, 
the outer and the jet regions. and are the additional wing 

vortex strengths on the wing surface outside and inside the jet, res- 
pectively. and arc the vortex strengths on the jet surface 

for the outer flow and the jet flow, respectively. In Eqs. (4.12) and 

(4.13) , 1‘i'Z and ^ I i z can be expressed as special forms of Eqs. 

(4.14) and (4.15). 
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If Eqs. (4 . 14)“(4 . 17) are introduced into Eqs. (4,10)-(4.13) 


these latter equations become 


On jet 
surface 


o - [ o i Cl t cc’JUiii - coi t,i = 

■^’4’^. , . u-n(i-yuO 


1 - ^ + 


U ■ e. 




i 

c s;;ji c ! + [ s ;;'3 u„ I - 1 ti- 4 [ so i c i - 

r 4 r 3 l«Cil = { - '^•+Tc/-'r:2^ \ 


On wing 
surface 


[CjiCi+rcji^./'i = 1 -^ 


±. ^ _ "2>4i 




U ► e. 


VC 

'•»Z 


\ 




c 3 { \ + c '^::'3 3 1 14 


^5 


I 




where the symbol ^ J represents a diagonal matrix. The system of 
Eqs. (4. 18)-(4. 21) can be solved either by eliminating some of the 
unknowns first and then inverting the resulting matrix equation, or 
by writing them as an augmented matrix equation for solution. The 
latter approach is used here to avoid great number of matrix opera- 
tions and great storage problems. For this purpose, Purcell's method 

A 

of solution [21] is used, because it requires only (Nj^^ + l)‘^/4 com- 
puter memory locations for the operation, where is the matrix size 
and much less computing time as compared with some standard direct 
matrix inversion methods. 


C^.2-0 



The resulting wing vortex distribution is assumed to be on the 
camber surface. Noting that the resulting pressure force is normal 
to the camber surface, the sectional lift and induced drag coefficients 
can then be obtained by resolving the pressure force in the proper 
direction and integrating over the local chord. ^ t sec- 

tional leading-edge thrust coefficient, then 

c 

^ . \ C u ♦ Cx) CC5 (o^ - Su>) Ax + c si V, (o< ~ 
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where 6 (x) is the angle between the tangent to the camber and the 
X-axis and same angle at the leading edge. Similarly, 

c 

(outside the jet) 

^ Cx) Ax “• (4. 

(inside the jet) 

In the above expressions, = ^w;u i ^*v/a. > where iCu is the wing vortex 
strength in the uniform flow. The sectional leading-edge thrust coef- 
ficient is computed by exactly the same manner as in a quasi vortex 
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lattice method in the thin wing theory (see Appendix C) . According to 
Eq.(47) of Appendix C, it is known that 


where C is the leading edge suction parameter. The expression for 
C can be derived from Eqs.(4.20) or (4.21) and similar equations for 
the uniform flow case. Aa an example, consider the portion of the 
wing in the outer flow. If the induced downwash at the i^^ leading- 
edge control point in the uniform flow is written as [N^]^ ^ 
then according to Eq.(46) of Appendix C, 

Wd, 7 = t \ i - \ 

Similarly, for the additional flow, it can be shown from Eq.(4.20) 
that 

5 U • t \ 7 

u-e^ ^ ^2. /i s 

where equal to - <x in the linear wing theory. 

With and C 2 computed, C is then equal to the sum of Cl and C^, 
because 

The sectional pitching moment coefficient about the y-axis can 
be found from Eqs,(4.27); 
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In Eqs.(4.27), the small contribution of the drag force to has 
been ignored* The over-all aerodyanmic characteristics of the wing 
are determined by spanwise integration of the sectional characteris- 
tics. For details, see Appendix C, 

To find the influence matrices in Eqs . (4. 14)-(4. 17) , it is assumed 
that the continuous wing and jet vortex sheets are replaced by step- 
wise-constant vortex distributions in the spanwise or the circumfer- 
ential direction. The induced velocity vector due to a "bounded” 
vortex element of strength Y dXg at any location has been found to 
be (see Appendix A) 
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where x measures the downstream distance on the vortex surface, and 
s ’ 
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The vortex element geometry Is shown in Fig, (4,1). Integration 
of Eq.(4.28) gives the total induced velocity: 



where the integration limits are chosen to cover the whole vortex 
sheet , 

The induced velocity due to the associated trailing vortices 
is, for one side of the strip, (See Appendix A ) 
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Other arguments in and have been omitted for simplicity. The 
subscript I denotes the spanwise or circumferential strip. The integral 


in Eq. (4.30) can be reduced to a finite sum by the quadrature method 
to be discussed below. (See Appendix C) 


to 





The integration interval is first transformed into (0, 1] through 


the relation 


at- 

where Xj^ is the lower limit and the upper limit. The results are 
then transformed into the 0-variable by 
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Hence , 
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This last equation is reduced to a finite sum by applying the mid- 
point modified trapezoidal rule [22]: 

i.(X,) = ^ iTce,) M, (e^, a:,) ^: ne^ 






where 


[_ I - cos &. ] 

6^ = ^ - c' =1. - ■ • . 

The integration points in Eq, (4.38) are specified by the modified 

trapezoidal rule, while the control points in Eq. (4.40) are dictated 

by the requirement of eliminating the Cauchy singularity contained 
1 

in Mjj . Similarly, Eq. (4.31) can be reduced to 

The above integration technique has been shown in Appendix C to 
work, well in thin wing theory where y is known to possess a square 
root singularity at x - . This is due to the fact that the square 
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root singularity has been cancelled by sin 6, so that there is no 
difficulty in applying the trapezoidal rule. When the integration 



interval is infinite, as is the case when the slipstream is assumed 
to be infinite, y is expected to decay in a certain way. Hence, a 
correct and efficient quadrature method to reduce the integration to 
a finite sum should be such that both the decaying nature and the 
Cauchy singularity in the integrand are accounted for simultaneously. 
Of course, such a scheme can be found by first removing the Cauchy 
singularity and then applying an appropriate quadrature method. 
However, it is expected that the consideration of infinite region is 
unnecessary in applications. Only that part of the jet surface which 
is close to the wing has to be considered. 

From Eqs.(4.37) and (4.41), the total induced velocity at i due 
to a vortex strip is 
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It follows that the normal velocity induced at i due to a unit vortex 
at k, or the (i,k) element of the influence matrix [N] , can be written 
as 
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where n^ is the unit normal vector at i. Similarly, 




The expressions for M, and are given in detail in Appendix D. It 
should be noted that if the influencing vortex is on the jet surface, 
8 should be replaced by 4 in Eqs.(4.43) and (4,44). 



5. Applications to Wing-Propeller Interactions 
5.1 Introduction 

As reviewed in the Introduction of Section 1, currently available 
methods in this area involve one or more assumptions in the formulation 
in addition to the linearization of the problem. For the purposes of 
presenting the objectives of the present applications, these assump- 
tions will be further discussed below. (1) Except Shollenberger , all 
existing methods used the Weissinger^s lifting line method or its 
variants for the wing representation. The distributions of sectional 
induced drag and pitching moment have not been obtained. (2) Ribner 
and Ellis [17] formulated the jet pressure condition in terms of velo- 
city potentials directly by integrating the pressure condition expres- 
sed in terms of tangential velocities. Further, they satisfied the 
pressure condition by using a doublet sheet associated with an unper- 
turbed flow field. Their results for the lift increment are comparable 
to those obtained by the slender body theory but are too low compared 
with Brenckmann^s experiments [23]. It is not known whether this is 
due to the inherent limitation of the linearized potential flow theory 
or due to the computational techniques. (3) No existing techniques 
are capable of accounting for different jet Mach number from the 
frees tream one, even in the linearized theory. Therefore, the high 
thrust condition where the jet Mach number may not be negligible 
can not be treated. 

With the def ficiencies of the existing techniques in mind, the 
objectives of the present formulation with applications to the wing- 



propeller interactions are as follows. (1) The linearized problem 
is to be solved as accurately as possible in order to determine whether 
the linearized potential flow theory is adequate for this type of 
problems. (2) Using the best available lifting surface technique, 
the and Cn, distributions can be determined and compared with 

available wind-tunnel data. In this way, if the linearized potential 
flow theory is to be used, the problem areas to be further studied 
can be identified. (3) Using the present lifting-surface technique, 
the importance of different Mach numbers inside and outside the slip- 
stream can be assessed. 

5.2 Mathematical Formulation 

The slipstream cross section will be assumed to be circular with 
constant diameter, and extended upstream as far as practical and down- 
stream to infinity, with its unperturbed velocity being constant and 
the same everywhere. Theoretically, the flow field ahead of the pro- 
peller due to the propeller alone can be approximated as that produced 
by a uniform sink distribution over the propeller disk [13]. Since 
this flow field is non uniform and is regarded as part of the base flow, 
not the perturbed flow, the perturbation due to the introduction of the 
wing is difficult to compute, because the governing partial differential 
equation would be one with variable coefficients. On the other hand, 
since the nonuniform flow inside the slipstream has been approximated 
by a uniform one which may be computed by the momentum theory, it is 
not well justified to use a more complicated model for the flow ahead 
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of the propeller. Therefore, the classical assumption of infinite 
silpstream is retained here. 

The right hand sides of the boundary conditions Eqs, (A.10)-(4.13) 
can be simplified for the following two cases: 

Case 1 . Slipstream not at an angle of attack. 

In practice, this may be the case where the propeller is not 
directly attached to the wing, a configuration which has been used 

^ A, — ^ A 

in some experiments. See Figure 5.1, Apparently, U • n = 0, 



Figure 5.1 Slipstream Geometry for Case 1 


u.l ^ V,. ; = V„ -hV, X. • = V, c. o< , where Is 

the rotational component of the propeller flow. It follows that Eqs. 


(4. 10)- (4. 13) can be written as 
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Note that the right hand side of Eq.(5,l) will be zero only if - M 


Case 2 . Slipstream at an angle of attack. 

This is the more realistic configuration where the propeller is 
fixed relative to the wing so that the slipstream will rotate with 
the wing. See Figure 5,2. In this case, it can be shown that 

A ^ ^ 

I' — i iMJsO -f- h ^ 

U=V.C</:SO(t 'f' ^iVi C< ^ 

/ •— * A i ^ ^ « 

yU := — yU «>< 




(i-, aa) 


(a~. ^ a ) 




Therefore, Eqs . (A,10)-(4. 13) become 
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Now, if the wing boundary condition in the uniform flow is taken to be: 
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Notice that Eqs.(5.5) and (5.11) have the same form if o< is small. 

Also notice the difference in wing boundary conditions inside the 
slipstream for the above two cases as defined by Eqs.(5.Aa) and (5.10a). 
The solution techniques described in Section A are used to solve the 
above problems. 

5.3 Numerical Results and Discussion 

Case 1 . Slipstream not at an angle of attack 

Some of the important existing methods are applicable only to 
this case. See, for example, Res. [14] and [17]. To check out the 
present program, a wing of A = 5.25 with a centered jet = 0.735 

as used by Rethorst [14] will be employed for comparison. No slip- 
stream rotation has been assumed. Instead of extending the slipstream 
to infinity both upstream and downstream, the present results were 
obtained by extending the slipstream 30 chord lengths upstream from 
the wing leading edge and 10 chord lengths downstream from the trail- 
ing edge for the purpose of satisfying the slipstream boundary condi- 
tions. This slipstream extent, was determined through numerical experi- 
mentation to ensure that the results will change as little as possible. 
Symmetry with respect to the x-z and x-y planes has been accounted 
for to reduce the number of unknowns. Only 8 circumferential vortex 





strips on the slipstream surface were used, as It was found that by 
increasing the circumferential vortex strips to 12, the results changed 
by only 0.5%. is assumed small so that cos c< 1 and sino( o< . 

The span loadings predicted by various methods are compared in Figure 
5.3. The results by other theories are taken directly from Ref. [17], 
except Shollenberger ^ s which are obtained from Ref [11b]. It is seen 
that the present method compares most favorably with Rethorst*s method 
[14]. It appears that both Ribner and Ellis, and Shollenberger pre- 
dicted lower loading outside the slipstream as compared with the 
present method. 

Figure 5.4 shows the comparison with the lifting line theory of 
Kleinstein and Liu for an elliptic wing of aspect ratio 6.5 with a 
centered jet. The nonuniform effects on the sectional lift curve 
slope have been included in this lifting line theory. Again, the slip- 
stream rotation is assumed zero in the computation. The lift increment 
inside the slipstream as predicted by Kleinstein and Liu is approxi- 
mately 10% higher than that predicted by the present method. 

Brenckmann's experimental results [23] for the span loading are 
compared in Figures 5. 5-5.6, where the sectional lift coefficient Cj^ 
is based on the average dynamic pressures. The experimental slip- 
stream rotation [23] was used in the present program by curve-fitting 
the measured results by a 4th-degree polynomial. However, the flow 
rotation outside the slipstream was assumed zero. Since the wing 
used in the measurement spanned the wind tunnel, in the present program 
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Figure 5.3 Comparison of Span Loadings Predicted by Several Theories. D/c=0.6 >^=0.735 
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an appropriate aspect ratio (A = 6.25) was chosen so that the program 
would predict correctly the sectional lift coefficient at the plane 
of symmetry for the wing alone case at c< = 6.5“. Again, 8 circum- 
ferential vortex strips were used on the slipstream surface. By 
increasing the circumferential strips to 10, the total lift coef- 
ficient was found to change by only 0.6%. The slipstream extent within 
which the slipstream boundary conditions are satisfied in the computa- 
tion was dictated by requiring less than 1% of change in the total 
lift coefficient for different slipstream lengths. In this way, the 
slipstream was extended 2.5 D upstream of the wing leading edge and 
5 D downstream of the trailing edge. In fact, the downstream extent 
was later found not as important as the upstream one. From the com- 
parison of the span loading, it can be seen that the agreement of the 
present results with the experiment is good in the region subject to 
the propeller upwash, but less satisfactory in the propeller downwash 
region. The total lift increment is compared in Figure 5.7, The pre- 
sent linear theory predicts the ratio AL/AL to be independent of oC , 
while the experiments show that it is a function of o< . The main reason 
may be as follows. At o< = 11.5®, experimental results indicated that 
some portion of the wing in the uniform flow was subject to flow 
separation. It follows that the lift coefficient in the uniform flow 
would be much lower than that predicted by the present potential flow 
theory. On the other hand, due to the blowing effects of the slip- 
stream, the potential flow results may be attained even at <x = 11.5®. 






Therefore, 4 L in the present method would be too low, because a too 
large Cj^ of the wing alone case has been subtracted. The results of 
Ribner and Ellis are seen to be closer to those obtained by the slender 
body theory. 

Case 2 . Slipstream at an angle of attack 

Even though this is the more realistic configuration in practice, 
no theoretical results for and distributions are available 

for comparison. It may be mentioned that Levinsky et al [18] have 
accounted for the slipstream inclination in their model. However, 
they satisfied the slipstream boundary conditions in the Trefftz plane 
and the wing tangency condition at the wing three-quarter chord line, 
as did Ribner and Ellis in Ref. [16]. Ribner and Ellis have found 
later that this procedure produced inaccurate results [17] . In addi- 
tion, no C 4 and c,„ distributions have been computed in Ref, [18]. 
Therefore, only comparison with experimental results will be made here. 

Recently, Nishimura of National Aeronautical Establishment of 
Canada [25] measured detailed aerodynamic characteristics of some 
half models with propeller diameter-wing chord ratio of A/ 3 and of 
NACA 4415 airfoil section. The propeller diameter-semispan ratio 
is 0,567 for the first wing to be compared. The measured pressure 
distribution was integrated to give the spanwise distribution of 

j C 4 3 nd Discrepancies exist in the integrated over-all charac- 
teristics Cj,and with the corrected balance measurements. The 

present comparison will be made with the pressure- integrated values. 



To use the present method for comparison with experiments, two 
difficulties are encountered. (1) The actual slipstream dynamic 
pressure distribution is far from uniform. Brenckmann defined an 
equivalent dynamic pressure for data reduction [23]. However, since 
no actual slipstream dynamic pressure distribution has been reported 
in Ref. [25] , an equivalent slipstream dymanic pressure can not be 
defined. In the following comparison, the slipstream dynamic pressure, 
and hence, the velocity ratio, is computed by the momentum theory as 
follows : 



where T^ is the thrust coefficient based on qg. (2) The slipstream 
rotational flow Vg was not reported. Levinsky et al [18] proposed a 
semi-empirical formula for Vg. However, this formula does not always 
produce correct shapes for the variation of the rotational angles, 
tan (Vg/Vg) , with radial distance as compared with some measurements, 
such as Ref. [23]. Since Brenckmann^s measurements of rotational angles 


are available, they are assumed to be applicable to other configurations 
if J is the same. Since ^ ^ > this criterion 

requires J C > to be the same. This criterion is obtained from 
the inviscid propeller theory with constant blade loading. For the 


present application with T^ = 0.8, J is found to be 0.316 [25]. From 
Eq. (5 .12) , yU can be computed as 0,44721. Therefore, = 0.5653. 

This same parameter is found to be 0.5667 for q* = 0.76 in Brenckmann's 
measurements. Therefore, Brenckmann's curve for tan~^ (Vg/Vg) at kt = 0.76 
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is used in the following computation. 

The camber slope for NACA 4415 airfoil can be described analyti- 
cally by the following relations [26] : 
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The slipstream extent to be included in the actual computation 
was again determined through numerical experimentation. In the fol- 
lowing numerical examples, the slipstream extent is taken to be from 
4 D ahead of the wing leading edge to 5 D behind the trailing edge. 
Since the experimental results are somewhat affected by the tunnel 
floor boundary layer and the nacelle, the following comparison is made 
only in terms of incremental changes. In addition, the center line 
of the slipstream is slightly below the wing plane in the experiment. 
This fact has been ignored in the computation. The interaction between 
slipstream surfaces have been included. However, the flow inside the 
slipstream on one side is not affected directly by its image region 
on the other side, but only Indirectly through the outer flow. Figure 
5.8 shows the spanwise distribution of ^ c and . The agree- 

ment is reasonable good in view of the assumption made regarding the 
slipstream rotation, a is seen to be too high on the propeller- 
downwash side. This is due to the under-prediction of the leading 
edge thrust. Some typical chordwise pressure distributions in this 
region and in the upwash region are compared in Figure 5.9, It is 
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Figure 5.8 Comparison of Predicted and Ac with Experiments for a 

Rectangular Wing of A=4. 705 at «=10’^ . D/c=4/3 , ==0. 44721 












seen that the predicted pressure level near the leading edge in the 
downwash region is too low. The negative loading near the leading edge 
is primarily due to the leading-edge camber. When the propeller down- 
wash was eliminated, the negative loading remained. However, with 
the propeller downwash included but the region x/ c iO.4 in the pro- 
peller downwash side replaced with uncambered section, the negative 
loading disappeared. This negative loading tends to shift the center 
of pressure backwards. The second feature of the theoretical 
distribution which may be mentioned is concerned with the values near 
any discontinuities in wing properties. The discontinuity in the 
present case is at the slipstream surface. Near the discontinuity 
outside the slipstream, becomes negative. This is due to the 
over-prediction of the leading edge thrust. Similar characteristics 
of the Ctj distribution occurs in wing theory in uniform flow for a 
variable sweep wing [27], The second wing compared is that of A = 3.22 
and full— span flap deflection of 15®. Most of the half span is covered 
by the slipstream with the propeller diameter-semispan ratio being 
0.83. The comparison of 4 and ^c^is presented in Figure 5.10, 
ACft^has not been compared, because the experimental distribution 
is difficult to read from the crowded curves. The incremental over- 
all aerodynamic characteristics are given in Table 5.1. 

In the applications related to Table 5.1, 245 simultaneous equa- 
tions have been solved with 34 K (decimal) memory and a total process- 
ing time of 11 minutes on Honeywell 635. No tapes have been used. 








Table 5,1 


Comparison of Predicted Aerodynamic Characteristics 
with Experiments with Wing-Slipstream Interaction 



A 

= 4.705 

A = 

3.22 


Present 

Experiment 

Present 

Experiment 

ACl 

0.771 

0.755 

1.414 

1.43 

ACd 

0.0752 

0.1 

0.3044 

0.39 


-0.1444 

-0.14 

-0.371 

-0.425 


Having established the accuracy of the present method in incom- 
pressible flow, it is of interest to see how the non-uniformity of 
Mach numbers would affect the predicted aerodynamic characteristics. 

For this purpose, the configuration used in Figure 5.8 will be used. 

The propeller thrust, and hence, the velocity ratio, is assumed constant 
while the freestream Mach number is increased from 0.1 to 0,3, The 
slipstream Mach number is taken to be proportional to the slipstream 
velocity. Let be the predicted sectional lift coefficient 

when the Mach number is assumed to be everywhere and the 

predicted sectional lift coefficient with Mach number non-uniformity. 

The percent increase, ( Mo>3 / C M*,) , is sho\im in Figure 

5.11. It should be noted that if the Prandtl-Glauert transformation 
is used locally, then 
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Therefore, for = 0.1, 0.2 and 0.3, Eq.(5.14) would give 2.1%, 

9.5% and 28.6%, respectively. It is seen that the sectional lift 
increment is much less than that predicted by the Prandtl-Glauert 
transformation. This fact can be explained as follows. As the slip- 
stream Mach number is increased, the sectional lift coefficient based 
on qg in the uniform flow with the corresponding Mg would also be 
increased. However, the effect of the slipstream surface, which is 
to reduce this uniform flow loading, is also increased. As discussed 
in Appendix B, this slipstream surface effect can be approximately 
explained by a reflection coefficient. In the two-dimensional flow 
this reflection coefficient can be written in the present notation as 
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If 3 ^— VQ/Vg = 0.44721, then ^21 = 0.667 if M^ = Mg. However, 

if Mq # Mg, then y\2l = 0.672, 0.691, 0.731, respectively, for = 0.1, 
0.2 and 0.3. Increase in the positive reflection coefficient implies 
increase in the reduction of loading inside the slipstream due to 
the slipstream surface, as discussed in Appendix B. On the other 
hand, the diffracted disturbances into the slipstream region from the 
outer flow, which are to oppose the reflection effects, are also slightly 
increased, with the values of the diffraction coefficient being 
0.7479, 0.7563 and 0.7741 for = 0.1, 0.2 and 0.3, respectively, 
as compared with 0.7454 for = Mg, where is defined as 
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However, because of the greater distance associated with wing vortices 
in the outer flow to the slipstream region, the diffracted disturb- 
ances in general have smaj-ler effects than the reflected ones, as 
discussed in Appendix B. Even though Eqs. (5.15) is strictly applicable 
only in the two-dimensional flow, it can be applied qualitatively 
to the present three-dimensional case. 

The percent increases in the total lift coefficient with a 
similar definition as given above are 0.52, 2,3 and 6.4, respectively, 
for Mq = 0.1, 0,2 and 0.3. Again, the increase in due to non uni- 
formity in Mach numbers is less than expected if the Prandtl-Glauert 
transformation is used, even though 56,7% of the wing span is immersed 
in the slipstream. 


5.4 Conclusions 

Within the assumptions of the linear potential flow theory, the 
wing-slipstream interaction problem cannot be solved by introducing 
a single slipstream vortex sheet in an equivalent uniform flow field 
if the slipstream Mach number is different from the freestream value. 
Instead, the perturbed flow fields must be solved separately for the 
slipstream and the outer flow regions and then matched on the slip- 
stream surface to satisfy the boundary conditions. For implementation 





of this concept in the computer, two vortex sheets on the slipstream 
surface have been introduced, one for the slipstream region and the 
other for the outer flow. This method also renders the Introduction 
of a source sheet unnecessary when the slipstream is at an angle of 
attack. Comparison with other analytical methods and experiments in 
incompressible flow showed that the present method is reasonably 
accurate in predicting the A and ^ distributions. In com- 

parison with experiments, some experimental distributions of slipstream 
rotational angles have been used in an approximate manner. It is 
recognized that for accurate prediction of aerodynamic characteristics, 
more accurate propeller rotational flow with both radial and axial 
variations is required, as the rotational flow changes the wing local 
angle of attack. However, for general applications, such precise 
characteristics of propeller rotational flow have yet to be established 
by either theoretical or empirical means. 

It was also found that nonuniformity in Mach numbers (with higher 
slipstream Mach number) increased the loading only slightly and to a 
magnitude much less than that which would have been predicted if the 
Prandtl-Glauert transformation was applied. 



6. Upper-Surface-Blowing Jet-Wing Interaction 
6,1 Introduction 

The current theory of a thin, jet-flapped wing was first formulated 
by Spence for a thin airfoil [41,42], Spence’s theory has since been 
extended to the three-dimensional wing with thin trailing -edge jet flaps. 
With the jet blowing over the flap, Spence assumed that the thin jet 
on the flap is "merely to prevent separation, and the pressure distri- 
bution is not modified in any other way" [42] . This assumption is reason- 
able if the jet has zero thickness. However, in the upper-surface blowing 
STOL configurations of current interest [7, 43], the jet sheet over the 
wing and flap is expected to be thick because of the use of high-bypass- 
ratio turbofan engines. The theory of thin trailing-edge jet flaps is 
not capable of predicting the high induced lift measured in some experi- 
ments (See "Discussions" below). . This is because the thick jet over the 
wing surface changes the indeed flow in much the same way as the solid 
ground would do to the wing flow. In other words, the jet which has 
higher dynamic pressure than the outer flow acts as a soft body which 
reflects the disturbances back to the wing to induce higher lift. 
Therefore, ignoring the thick jet effect becomes unrealistic in this 
application. 

In the following, a linear theory for a thick, jet-flapped wing will 
be formulated, following the general theory of wing-jet interaction 
presented in Section 4. With this theory, the aerodynamic characteristics 
of wings of arbitrary aspect ratios with jet interaction can be predicted. 



The effects of jet thickness, jet chordwise, spanwise or vertical location, 
Mach number nonuniformity and jet temperature have been included. The 
effects of jet deflection are accounted for in a similar manner as has 
been done by Spence [41] . The basic assumptions are (1) linear and 
inviscid theory, (2) no jet entrainment, (3) constant jet thickness and 
velocity and (4) no wing thickness and nacelle. 


6,2 Mathematical Formulation 

In the present formulation, the jet is assumed to be of rectangular 
shape located at an arbitrary vertical distance from the wing surface. 
Therefore, the wing is assumed to be completely in the outer flow even 
if the jet lower surface may be on the wing upper surface. According to 
Eqs , (4 , 10) -(4 ,12) , the boundary conditions to be satisfied for the 
additional vortex strengths are 
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Again, using the vortex sheets for the jet surface, Eqs. (6, l)-(6.3) can 
be written in terms of the influence coefficient matrices as follows: 



On jet 
surface 


(6.4) 


^_U- n (l-y^O -^ 4 . 

* u-£. 




1 - 




S 


i- T.^7 


^3. 


] 


(6.5) 


(on wing surface) 

Eqs. (6.4)-(6.6) are written in the form of an augmented matrix 
equation and solved in the manner as discussed in Section 4. Some special 
features of techniques used in solving Eqs. (6.4)-(6. 6) are discussed 
below. 


(1) That part of the jet which is directly on the wing surface is 
assumed to be subjected only to the free stream angle of attack, dis- 
regarding whether there is flap deflection or not. The jet deflection 
due to Coanda effect will be discussed later and does not change U » n 
term in Eq.(6.4). Therefore, the first term on the right hand side of 
Eq.(6,4) can be written as, for the jet control points not on the sidewalls. 


Ll ( i - mO r ^ t 

= - L t U-y^')t««o<l 


(6.7) 


where the upper sign is for the upper jet surface. On the other hand, 
for the jet sheet behind the trailing edge when the jet is deflected 
because of flap deflection, the jet will be subjected to varying angle 
of attack due to the continuous changing of the jet sheet curvature. 


C, i 



This curvature distribution is unknown and must be solved simultaneously 
with the vortex strengths. To avoid Iteration, Spencers method [41] 
can be adopted here in the following way. Since the jet flow is irrota- 
tional, it is true that (See Figure 6,1) 


M(R- V, (Rt it,) 

or 


(6.8a) 


V, - V = 

where Vj = (vj^ + V2)/2. Note that both vj^ and v 
jet velocities. Introducing the nondimensional 
as given by Eq.(4.2), it is seen that 


(6.8b) 

2 are dimensional total 
perturbed velocity potential 



Figure 1. Geometry Illustrating Irrotationality of a Deformed Jet 



since Eq.(6.2) can be written as 
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by applying Eq.(6,2a) to the upper and loitrer surfaces of the jet element, 
it can be shown that Eq,(6.9) is reduced to 


^ “ 1 ^ ( 6 . 10 ) 
Eq. (10) is seen to be consistent with Spence’s results when the thin 
jet assumption is made. In other words, if then 

* J~>- i 

the jet vortex strength would be =r > consistent with 

Eq.(26) of Ref. (41), Assuming that ^ ^ j as has been 

done in Eq.(4.8), Eq. (6.10) can be reduced to 










(6.11) 


For convenience, the right hand side of Eq,(6,ll) will be denoted by f(x): 


( 6 . 12 ) 


= T - (^X+ 

For small jet deflection, 1/R = d^z/dx^ . It follows that 

X ft. = 

■j|Ct.6.) = (6.13) 

where 6^ is the flap angle. Eq. (6.13) can be integrated to find the jet 
slope ^ at an arbitrary point. Once ^ is obtained, n becomes 



where the upper sign is for the upper surface. Notice that f(x) in Eq.(6.12), 
dz 

and hence — , depends on the unknown jet vortices These terms 

can be incorporated to the left hand side of Eq. (6,4) before the solution 
is attempted. 

To integrate Eq.(6.13), two methods have been used. The first is 
the point-slope formula which is the following: 


(6.16) 

The second method is as follows: Let 

where Cj is the jet length from the wing trailing edge to be included 
in the analysis. Then 


or, 






It follows that 


o ^ 


(6.18) 



where the conventional trapezoidal rule has been used. in Eq.(6.18) 

are seen to be exactly the angles used in defining the jet control points. 
The second method has been found to be more efficient than the first one. 


(2) With the jet on the wing surface, the evaluation of the u-induced 
velocity, and therefore [S] matrices and in Eq.(6.5), 

must be handled carefully. To illustrate this point, consider the expres- 
sion for u(x,z) for the 2-D case in incompressible flow: 


U(V 2)= — (6-19) 

Observe that for 0<x<l the integrand in Eq. (6,19) has a second order 
singularity as z approaches zero while the whole integral will be multiplied 
by z. As has been checked numerically, Eq. (6.19) can be integrated accur- 
ately by the integration technique discussed in Appendix C, if z and the 
number of integration points are not too small. Since in the present 
integration technique, just like the conventional VLM, the number of 
integration points is equal to the number of unknowns to be solved, 
increasing the number of integration points becomes numerically unrealistic. 
On the other hand, as has been illustrated in Appendix C, the present 
technique is quite efficient in treating integrals with Cauchy singularity. 
Therefore, the best numerical technique under the present circumstances 
is to rewrite Eq.(6,19) as 
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The last integral in Eq,(6.20) can be integrated in a close form: 
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Eqs.(6.21) and (6.20) show that u(x,z) becomes y(x)/2 if z = 0. However, 
in the three-dimensional case, the corresponding integral is too tedious, 
if not impossible, to Integrate. Therefore, the following numerical 
method of integration has been selected. The integral in Eq, (6.21) 
can be integrated accurately by the method illustrated in Appendix C 
by increasing the number of integration points without changing the 
number of control points in such a way that the integration points are 
always interdigitated with the control points (so that they never coincide) . 
Intuitively the number of integration points N can not be arbitrary, 
instead should be a certain multiple of the number of control points M. 

The specific relation between the tv7o can be derived as follows. As 
shown in Appendix C, 




To (6.22) 

where 0]^'s are roots of Tj^(cos 0j^) = cos N0j^ = 0. If 6 is chosen so that 
Tjq*(cos 0) = 0, then 6 and 0j^ will never be the same, because the zeros 
of Tjj(cos 0) and Tjj"(cos 0) never coincide. Since Tj^'*(cos 0) = N sin N0/sin 0 
(See Appendix C), it follows that 0 should be such that 

sin N0 = 0 (6.23) 

By trigonemetric relations, Eq.(6.23) can be written as 



N0 


Thus, Eq.(6,23) can be satisfied if sin 0 where p is any integer 

It follows that 0 for interdigitat ion should be such that 
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Let N/2^ = M, the number of control points. Therefore, the control 
points are given by 0^^ = iw/M, the same as has been found in Appendix C, 
With M as the number of control points, the number of integration points 
(N) for interdigitat ion is, therefore, 

N = 2% (6.26) 

and 
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As has been checked numerically, the error in numerically integrating 
the integral in Eq,(6.21) is less than 1% for p = 3 (i.e., N = 8M) , 
z = 0.05 and 2 ^ M ^ 5. 

For the three-dimensional case, the above technique is used only 
for the integration strip with the control point directly below or above 
it. All the other integration strips are reduced to a finite sum by 
the method described in Section A, For the particular strip under con- 
sideration, the u-induced velocity, Au, due to this strip of vortex 
distribution is written as 
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The expression for • x Is given in Appendix D. For unswept bound 
vortex distribution, the integral ) M • c can be easily integrated 

in a closed form. This provides another check case for the accuracy 
of the present numerical technique. This is done in Appendix E. It 
may be mentioned that the increase in computing time due to the refine- 
ment in computing the u-induced velocity is minimal, but the increase 
in accuracy is considerable, 

(3) When the jet is on the wing surface, some of the jet control 
points coincide with the wing control points. At these particular jet 
control points, the sum of the first two terms of Eq.(6,4) would be 
zero because of Eq.(6,6), i,e,, 

Therefore, if the system of equations is solved by regarding 

wa 

{Yoj^f unknown vector in that order and process Eq,(6.6) 

first and then Eqs. (6.4) and (6.5), as has been done in Section 5 for 
wing-slipstream interaction, difficulty would occur in the solution 
method. As a simple example, consider the following system of equations: 

X + 2y = 0 (6.29a) 

X + 2y + z = 1 (6,29b) 

X + y + z = 2 (6,29c) 

If X and y are to be solved first from Eqs, (6.29a) and (6.29b) in terms 

of 2 to be substituted into Eq. (6,29c), great difficulty would be encountered. 

On the other hand, if z is regarded as the first component of the unknoxim 



solution vector and solve it first, the difficulty would disappear. 

Using this concept, the unknowns in Eqs , (6, 4)-(6. 6) are now solved in 
the following order: ^^wa^ * Eq-(6.5) is processed first 

and then Eqs. (6,4) and (6.6), It appears tliat this order is also natural 
in that the u-induced velocity due to unit vortex strength for the jet 
behind the trailing edge is required in computing dz/dx for the trailing 
jet when the jet deflection is present (See Eq,(6,ll)), This information 
can be stored on tape to be used in Eq.(6.4). 

6,3 Numerical Results and Discussions 

Since no applicable experimental results are available for comparing 
detailed aerodynamic characteristics, the computer program was slightly 
modified to treat the conventional jet flap with thin and stiff jet 
assumptions. By stiff jet, it is meant that Vj/Vq -J- » but Vj^t^ is finite, 
where tj is the jet thickness. It follows that p = 0 and the jet flow 
is unperturbable. This implies that 3q'^/9n = 0 in Eq.(6.1). Furthermore, 
Eq.(6.10) gives 

R ^ ^ (6 30) 

where cp (y) is the sectional jet momentum coefficient. To avoid difficulty 
with the presence of the jet sidewall vortices, the jet thickness was 
assumed to be 5% of the chord length. Das’ theory and experiments [44] 
were chosen for comparison. To investigate the numerical convergence, 
various lengths of the trailing jet sheet were used until a maximum value 


6.II 







for the predicted Cj^ was reached for a given number of streamwise vortices. 
It was found that the jet sheet curvature, and therefore the jet vortex 
strength, decreases rapidly downstream from the wing trailing edge. The 
results for a wing with half semispan blowing with c^(y) = 1.0 and 30“ 
jet deflection are shown in Figure 6.2, The results indicate the cor- 
rectness of the present computer program and the mathematical formulation 
for the trailing-edge jet flap. 

As the two-dimensional computer program used in Appendix B is 
believed to be correct, it is used again here to check the numerical order 
of magnitudes of the predicted sectional lift coefficient by the present 
3-D program. For this purpose, Phelps^ wing planform is used, except 
that the leading edge Krueger flap and the trailing edge flap extension 
are ignored. This configuration is shoxTn in Figure 6.3. The camber 



Figure 6,3. Geometry of Wing Planform for Numerical Investigation 



and the jet angle of attack are set to zero for the comparison with the 
2-D theory. At the spanwise station 2y/b = 0.26667, the jet thickness 
ratio tj/c is 0.1147. For a velocity ratio Vj/V^ = 5.7314, the 3-D results 
show higher lift augmentation (See Table 6. iXphysically, this may be due 

to the fact that the disturbances diffracted into the jet region have 
been weakened because of 3-D spreading effect and therefore they become 
weaker as they are diffracted back to affect the wing flow. According 
to the 2-D theory (Appendix B), the diffracted disturbances are to 
decrease the wing lift and therefore, weaker diffracted disturbances 
would increase the lift. 


Table 6.1. 

Comparison of 2 

-D and 3-D Jet 

Lift Augmentation. o( = 10"^ 

, t./c = 0.1147 

’ J 


No Jet 

With Jet 

3-D 

0.78876 

4.7465 

2-D 

1.0966 

c 4.487 

18.5598 (No diffractions) 


The experimental results of Ref. 7 showed only the total aerodynamic 
characteristics. To compare with them, the jet velocity would be needed 
in the present method. This is computed by the momentum theory as 
follows. The thrust is given by 




(6.31) 



Then, 
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Solving for V./V , it is found that 
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For the jet thickness of 2 inches (0.1667 ft), Aj = 0.2167 sq. ft. 

The half wing area used by Phelps, et al, in computing the thrust coef- 
ficient C was 5.61 sq. ft. Using this information, Eq.(6.32) gives 
Vj/Vq = 5.7314 for = 2.095 and pj = . 

Another difficulty must be solved before the comparison with experi- 
mental results can be made. As the flap is deflected, the jet will 
follow the flap surface due to the Coanda effect. The Coanda effect is 
due to the high suction on the flap surface to attract the jet downward. 

It exists also under wind-off conditions. Since this effect can be explained 
only by the viscous theory, its exact prediction is out of scope of the 
present investigation. For the present purpose, it will be estimated 
in accordance with the momentum theory. As illustrated in Figure 6.4, 
when the jet is deflected at the flap knee, a jet reaction F^ equal to 

H 5- Af V*" /-) produced. Its 

lift and drag components are therefore 
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(6.33) 
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(6,34) 
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In coefficient form, they become 


= 1- i +c<) (6.35) 

Cj,_^ = ^ s;»i(S(/ 2 ) sin(S ,:/2 t •=<) (6.36) 

Under wind-off conditions, Eq. (6.33) gives = Tj- sin 6^ which is 
consistent with the equation given Glahn [45]. However, the expression 
for Dj, suggested by Glahn is Dj. = sin $£ tan 6^ , while Eq. (6.34) 
would give Dj, = 2T sin^ 0^) , 



Figure 6.4, Coanda Effect 



For = 30°, a = 5° and = p j , Eq, (6,35) gives the Coanda lift 
coefficient as 1.2277. The aerodynamically induced lift increment (ACl)jI^ 
which is the difference between the augmented lift coefficient without 
the Coanda effect and the wing alone case, is plotted in Figure 6.5 
with respect to the jet leading edge location. The jet exit in Ref. 7 
is estimated to be at 1/3 of_the wing chord (excluding the flap) from 
the leading edge. That is, = 0,7 x 1/3 = 0.233, From 

Figure 6.5, (ACl)^^ is found to be 1.8. Therefore, the total lift incre- 
ment is ACl “ 3.0277(=1.8 + 1.2277). This compares favorably with 2.9 
from Figure 12 of Ref. 7. From Figure 6.5, it is also seen that the 
trailing-edge jet alone will not produce enough lift to match the measured 
value. Since the present method gives only the induced drag, while 
the experimental results included all drag components, the drag computa- 
tion is not compared. The experimental pitching moment is directly 
affected by the thrust effect, because the engines are above the moment 
center. Besides, the Krueger flap and the flap extension may also affect 
the change in the pitching moment. Therefore, the pitching moment has 
not been compared. In Figure 6.5, the predicted span loading has also 
been shown. The incremental loading is seen to be concentrated in the 
jet region. The higher near the jet boundary is probably due to the 
jet sidewall effect which allows the diffracted disturbances in the jet 
region to leak out through the sidewall so as to reduce their effect on 
the wing flow and also due to the increase in the jet thickness/chord 
ratio outboard. In reality, the loading will probably fair more smoothly 






into the outside region because of the jet sideward spreading effect 
and the turbulent mixing. 

Putting = 0, additional data are generated for the configuration 
of Figure 6.3 for the effects of the jet vertical location on the 
aerodynamic characteristics. Table 6.2 shows the change in the total 
aerodynamic characteristics at three heights. The rapid change in Cl 
and Cq] with jet height seems to confirm the 2-D results of Appendix B. 


Table 6.2. Change in Aerodynamic Characteristics 



due to 

Vertical Shift 

of Jet Location. 



o 

> 

> 

= 5 .7314, o( = 

5% Mo = Mj = 0 


h, ft . 

h/tj 



AX 

r 

0 

0 

0.7484 

0.04715 

-0.1417 

0.07 

0.42 

0.1545 

0.01357 

-0.07298 

0.15 

0.90 

0.0487 

0.0057 

-0.01156 


The appropriate parameter here is the height /thickness ratio, or h/tj . 
Large velocity ratio tends to increase the rate of decrease in 

ACl with respect to k/tj. Also, the aerodynamic jet interaction is to 
shift the center of pressure forward. Note that this does not include 
the effect of the Coanda jet reaction when h = 0. This shift in the 
center of pressure is also evidenced in the pressure plot of ’Figure 6.6. 
The span loading and the induced drag distribution are compared in 







Figure 6.7, Figure 6,8 shows the effects of mach number nonuniformity 
on the span loading. The jet mach number is computed by the relation; 



and T = p^/pj. The general trend is similar to that predicted by the 
2-D theory. That is, higher jet mach number produces higher lift and 
high jet temperature tend to degrade the lift capability. 

As a final item for discussion, it is of interest to note that Eq,(6.1) 
shows that when the jet is subject to an angle of attack, a source dis- 

A 

tribution would be required on the jet lower surface where n = -k; while 

on the upper surface a sink distribution would be needed. Since the 

jet lower surface is closer to the wing, the net result is a reduction 

in lift because the source distribution produces a net downwash on the 

wing. This observation has been checked out numerically on the computer. 

This interaction characteristics may have the following effect. As the 

wing is subject to a positive roll rate, the right wing will see an 

increase in angle of attack, thus increasing the lift. However, the 

increase in angle of attack for the jet will reduce the lift as discussed 

above. The net effect must be checked to see how the roll damping 

derivative will be affected. Similar remarks are also applicable to 
P 

the pitch damping derivative This effect has not been investigated 

in the present study. 

In the above numerical investigation, a typical computer run required 
33 K (decimal) of memory, 5.5 minutes of processing time (6 min. with Hach 
number nonuniformity) with Honeywell 635 computer solving 198 unknowns. 

One tape is needed. 
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APPENDIX A 


Derivation of Induced Velocities Due To A 
General Vortex Segment in Linearized Compressible Flow 


The velocity field due to a line segment of vortex of strength V 
is (p. 43 of Ref. 28 ) 
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Substituting Eq. (5) into Eq.(l), it is found that 
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Eq, (15) can be recast to a different form without involving ^ and £/ 

■— -J * - ^ _. ^ 

Since ,£ = b - a , and £' = b ^ - a ' , it can be shown that 
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Eq, (15) becomes 

In incompressible flow, ^8 = 1. In this special case, Eq. (16) becomes 
the expression used by Maskew [ 29 ]. 

For the induced flow field due to the trailing vortices, Eq. (1) 
must be used, where L now starts from the bound vortex position Xg 
to CO , Thus, the induced velocity due to one trailing vortex on one 
side of a bound vortex is given by 
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APPENDIX B 


Compressibility Effects on Airfoil-Jet Interaction 


1. Mathematical Formulation 

Since the general three-dimensional problem is much more dif- 
ficult to analyze for the Mach number effect, a two-dimensional an- 
alysis for airfoil-jet interaction is first undertaken. A plane jet 
with Mach number M^ is assumed to be imbedded in a freest ream of Mach 
number Mj^. A thin airfoil is placed at a distance Z' below the lower 
jet surface, as shown in Figure 1. This may represent an idealized 
configuration of upper-surface blown jet flapV'"*X" similar configura- 
tion has been analyzed in Ref. [11a] for incompressible flow. The 
case where the airfoil is in the mid-stream with different velocity 
regions above and below in incompressible flow has been treated by 
Ting and Liu [30]. Since the thin airfoil in the stream may be re- 
presented by a vortex distribution, the expression for the induced 
velocity vector due to a two-dimensional vortex in the linearized 



Figure 1. 


Geometry Definition 



compressible flow must be derived. Noting that the vortex axis is 


in the y-direction, Z ~ (yi-y, Eq,(4.28). From Eqs.(4.29) 

of Section 4, the following relations may be found: 
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It follows that if ?TdXg =. p , Eq. (4.28) gives, for a two-dimensional 
vortex, 



Hjir 

‘fTr 


i ( ^ - zO C - I ^ ^ - y, ) i. 


/z - y'X 

[ (a- ^0 V^''( Xx- y)+ -ii ^ ^ 


P^(y, ~y)(y.-y.) ) 

[ f ;jc - at 0^ + " ( X. - y /i- /S " ( 2 - Z o ^‘3 i 

d Cz- z9 - (X-xO ^ ^ P(3 

(X-X'y+p-^{Z-Z‘T ^ <-X-X'T t p\z-Z‘T 


To examine how disturbances reflect and diffract at the inter- 
face of two flows of different energy levels, consider a vortex of 
at z ' = -h in region 1. See Fig. 2. The boundary conditions at 
z ~ 0 are 
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(Pressure continuity) 



Due to the presence of the vortex at (x* , -h) , the flow 
field In region 2 will be disturbed. According to Koning's concept 
[12] for incompressible flow, the disturbed flow field in region 2 
due to the vortex can be represented by an additional vortex which 
must be situated outside the region. In the present case, a vortex 
of strength b at z' = be introduced. It follows that the 

induced flow in region 2 can be written as 



Figure 2. Vortex Images 
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On the other hand, the disturbance due to the vortex will be re- 
flected back to region 1 at the interface. This reflected disturbance 
can be represented by a vortex of strength "a" at the image point 
in region 2. Hence, the total induced flow in region 1 is 

_ ll^i e(zfh)- ?;(>.-;>:< ) , A ^ ij -A.') 

Using Eqs. (4) and (5), the boundary condition Eq. (2) becomes 
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Similarly, the pressure condition becomes 
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Solving Eqs. (6) and (7) simultaneously for a and b, it is found that 
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Eqs, (8) and (9) show that the parameters determining the additional 
vortex strength *'a” of the reflected disturbances and the vortex strength 
"b" for the diffracted disturbances are the dynamic pressure ratio, 
the velocity ratio and the ^ ratio. It is seen from Eq. (6) that the 
ratio of the vortex strengths in the two regions depends not only 
on the velocity ratio but also on the ^ ratio. In incompressible 
flow, Eqs. (8) and (9) become 
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which can be seen to be the same as those given by Koning [p. 390, 
Ref. 12]. 


The above results can be summarized as follows: (1) Reflection 

is accounted for by putting an additional vortex with strength "a" 
at a distance equal to that of the originating vortex from the inter- 
face and at the opposite side of the interface. (2) Diffraction is 
accounted for by putting a vortex of strength ”b” at a distance equal 
to that of the originating vortex from the interface multiplied by 

'®(incident)^ /^(diffracted) “rginating 

vortex. 

Now apply the above principles to the configuration shown in 
Fig. 1. Define and Z 2 as the distances from the lower and upper 
jet boundaries, and B 2 , respectively. The distance will be regarded 
as positive if the vortex is above the boundary. Let the "reflection" 
and "diffraction" coefficients A. and S' be defined as follows: 
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A vortex / at a distance h from (or Zj^ = -h) is reflected at 
with an additional vortex A.,^ ^ at Z-^ = +h and diffracted into the 
midstream with a vortex at h. The diffracted distur- 

r'Jjrlec.+Mi.J o 

bance is then^at B 2 to give an additional vortex ^ ^2 “ p h + t. 

The reflected disturbances are in turn reflected and diffracted at 
Bj^. The diffracted disturbances into the lower region are due to the 
vortex ^ at Z, == , and so on. 

These are illustrated in Fig. 3. If (^\ ~ (*>i, » the first diffracted 
image for the lower region would be at 'Z^— with strength 

equal to 




( V.-V,"., 

^ 2-V, V. \ 




This last expression agrees with that given in Eq. (10) of Ref. 30. 
(Eq. (10) of Ref, 30 seems in error in the image locations, as can 
be easily checked. The originating vortex is in the upper stream 
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while the images are also in the upper stream. This violates the 
principle that no additional flow singularities can be introduced into 
the region under consideration [12]), Adding the effects of all image 
vortices in the upper stream gives the total diffracted disturbances 
in the lower stream. With reference to the airfoil coordinate system, 
the downwash produced at the airfoil is given by 
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where the downwash is normalized with respect to and all lengths 
are based on the chord length. The airfoil integral equation be- 
comes, using Eq.(16), 
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Following the classical airfoil theory, Eq.(17) can be solved 
by assuming 
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when Eq. (18) is substituted into Eq, (16) , some integrals must be 
evaluated. It is known that 
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From Ref. (30), it is obtained that 
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the airfoil integral equation can be reduced to 
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Eq,(30) can be solved by expanding all ’s and <l^c/dxin cosine 

Fourier series and equating the coefficients of cosine terms on both 
sides. For this purpose, the program listed in Ref. (31) was revised 
to solve Eq,(30). 


2. Numerical Results and Discussions 

To check the correctness of the computer program, Shollenberger^ s 
results [11a] for incompressible flow were compared. According to 
Figure 10 of Ref . [11a] , for an actuator disk thrust coefficient = 0.5 
with, the airfoil just at the lower jet boundary, is found to be 
1.44 at o< = 10®. This is the lift coefficient for the flat-plate 
airfoil and simulated nacelle combination with the airfoil one chord 
length behind the nacelle. In the present program, no nacelle was 
included and the jet was extended to infinite upstream. Furthermore, 
no jet boundary displacement was assumed. To use the program, the 
velocity ratio was needed. This was computed by the momentum theory 
as follows : 
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The program gave <^j^/a7to< = 1.2016 at Vy ^V, ) ” 1.225. There- 

fore, may be computed to be 1.3177. This result is seen to have 
the same order of magnitude as that given in Ref. [11a]. 

The results for some incompressible cases for t/c = 0.5 are 
presented in Fig. 4. The variation of the lift and moment coefficients 
with the velocity ratio and the vertical distance is typical. As 
the airfoil is moved away from the jet, both and are de- 

creased. On the other hand, if the velocity ratio is increased, both 
and — are also increased. The dashed lines represent the 

case where all diffracted disturbances are ignored, so that only the 
first two terms in Eq, (16) are included in the analysis. The results 
show that the reflection represents the primary contribution in the 
jet interaction for this jet thickness ratio. Notice also the center 
of pressure shift due to the interaction. The interaction tends to 
shift the c.p. backward until the airfoil is moved close to the jet 
surface. By comparing the two curves in Figure 4c, it may be concluded 
that diffracted disturbances tend to move the c.p. forwards and re- 
flected disturbances would move the c.p. in the opposite direction. 
Since the diffraction terms are more important as the airfoil gets 
closer to the jet boundary, this explains why the c.p. will shift 
rapidly forwards as the airfoil is moved toward the jet boundary. 
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Figure A Variation of Aerodynamic Characteristics of a Flat- 
Platc Airfoil with Jet Interaction in Incompressible 
Flow. t/c=0.5 
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The values of the "reflection" and "diffraction" coefficients as 
defined in Eqs. (12)-(15) are given in Table 1. Notice that the in- 
crease in c^ and— is associated with more negative "reflection" 
coefficient x 

In the compressible flow, T 2 /Tj^ = 3.0 and = 0.4 were assumed. 
The Mach number ratio and the density ratio are computed by the fol- 
lowing relations: 
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Table 1 

Reflection and Diffraction Coefficients 

Case 1 Case 2 

= 0, T 2 /T 1 - 1.0 = 0.4, T 2 /T 1 = 3.0 


V2 / V1 

3.0 

1.5 

3.0 

h 2 

- 0.8 

0.1542 

- 0.5845 

^21 

0.8 

- 0.1542 

0.5845 

"^12 

0.6 

1.6915 

1.5845 

hi 

0.6 

0.5771 

0.4155 





Eq, (33) implies that pj^ = P 2 » The results are shown in Figure 5, 

The diffraction is seen to be important for ~ unimpor- 

tant for V 2 /Vj^ = 1.5. Figure 6 shows how the nonuniformity in Mach 
numbers will affect the aerodynamic characteristics. It is seen 
that the loading is slightly increased for the higher jet Mach number 
with other parameters remaining unchanged. This fact can be explained 
by the more negative reflection coefficient. For T 2 /T-|^ = 1.0, 

V 2 /V 1 = 1.5 and Mj^ = M 2 = 0.4, computed to be -0.3846. 

On the other hand, if M 2 = Mj^(V 2 /V 2 ) = 0.6, then A -22 “ -0.4410, 

More negative reflection coefficient implies more lift. Since 
remains to be -^0.3846 even if = M 2 = 0, the above results indicate 
that the lift augmentation, expressed as compressible 

flow is higher than that in incompressible flow. 

Cross examination of Figures 4-6 reveals the following results. 
Comparing Figures 5 and 6 for V 2 /Vj^ = 1.5, it may be concluded that 
high jet temperature has detrimental effects on loading and produces 
less nose-down pitching moment about the leading edge. At ~ 

and V 2 jVj = 1.5 , is less than the uniform flow value because ^12 
is now positive. Examining Figures 4 and 5 for ~ 3.0, it is 

also seen that the lift augmentation for T 2 /Tj^ =3.0 and = 0.4 is 
much less than that in the incompressible flow when the airfoil is 
close to the jet boundary. Therefore, wind-tunnel tests of wing-jet 
interaction with cold jet may overestimate the loading unless the 
dynamic pressures are correctly simulated. 
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For simplified quantitative estimation of the wing-jet interaction, 
it has been mentioned above that the reflection accounts for most 
of the contribution under the specified conditions. In fact, a criti- 
cal thickness ratio may be defined as that where omission of the dif- 
fraction terms would incur less than 1 % error. This critical thick- 
ness ratio depends on V2/Vj^, T2/T1 and Two typical curves are 

plotted in Figure 7 for airfoils situated at the lower jet boundary. 

It is seen that for a hot jet in a compressible external flow, the 
interaction effect may be estimated quantitatively by accounting only 
for the reflected disturbances without incurring too much error. Of 
course, in general, the reflection coefficient is still qualitatively 
a good measure of how the loading will change with the change in some 
of the physical parameters. 

One interesting result applicable to the upper-surface blowing 
concept may be mentioned. This is based on the following observations. 
( 1 ) With negative ^12 5 i.e., the lift is increased over the uniform 
flow value, diffraction terms in Eq. ( 16 ) are to decrease the loading. 
These diffraction terms become unimportant as the jet thickness is 
increased. ( 2 ) The loading is always increased when increased. 

Keeping these two facts in mind, now consider the case with a given 
amount of thrust. As the jet thickness is increased, the diffraction 
becomes less important so that the lift is increased, while the de- 
crease in velocity ratio is to decrease the lift. On the other hand, 
as the jet thickness is decreased, the above two effects reverse the 
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trend. However, for h “ 0 and t = 0, the diffraction terms are so 
important that they exactly cancel the reflection effect so that the 
uniform— flow results are obtained. This can be seen by summing the 
coefficients inside the curly brackets in Eq.(l6). This gives 
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But for incompressible case, 
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It follows that 

f = 1 t ^|X ^2t ^ , 

independent of the velocity ratio. This result is also true in the 


compressible flow. Therefore, as the jet thickness is increased from 



zero, the lift will be increased. Eventually, the lift will be de- 
creased after reaching a maximum value because of the low velocity 
ratio associated with the large jet thickness. An optimum jet thick- 
ness is seen to exist. For the purpose of numerically computing the 
optimum jet thickness, a thrust coefficient of 2.095 per unit 
span is assumed, a value taken directly from Ref. [7]. This thrust 
coefficient is defined as the amount of thrust divided by the jet 
thickness and the freestream dynamic pressure. By the momentum theory, 
the velocity ratio then depends on the jet thickness in accordance with 
the relation: 
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The results of computation are presented in Figure 8. It is seen 
that an optimum jet thickness does exist. For T2^"^l ~ 1*0, the op- 
timum thickness ratio t/c lies between 0.175 and 0.2, approximately 
independent of the Mach number. On the other hand, for ’^ 2 ^’^!. ~ 3.0, 
the optimum t/c becomes 0.1, Some assumptions made in deriving this 
result are as follows. (1) The jet has been extended upstream to 
infinity, even though in reality the jet may be issued from the air- 
foil surface. This fact may affect the results quantitatively. How- 
ever, the qualitative conclusion is expected to remain unchanged. 

(2) The jet static pressure is assumed to match the ambient static 
pressure. (3) No turbulent mixing effect is included. (A) No air- 
foil thickness effect has been accounted for in tVie present linear 
theory. 




Even though the present numerical results and discussions are 
related to one specific configuration only, the following qualita- 
tive generalization is valid. Positive reflection coefficient implies 
decrease in lift from the uniform flow value. The contrary is true 
for a negative reflection coefficient. Reflection of airfoil distur- 
bances generated in the high dynamic pressure region from the jet 
surface into the same region is to decrease the airfoil lift from the 
uniform flow value. This case is signified by a positive reflection 
coefficient. For example, if the airfoil in Figure 1 is placed in 
the jet, the uniform flow conditions are then the jet conditions. 

The reflection from the jet interface will decrease this uniform-flow 
loading. On the other hand, reflection of disturbances generated in 
the low dynamic pressure region from the jet surface into the same 
region is to increase the loading from the uniform flow value. This 
last statement applies to the configuration of Figure 1. In general, 
the diffraction effect is to oppose the reflection one. However, 
in most applications, the diffraction effect is small. The above 
generalization is very useful in making interpretation of three-dimen- 
sional results in wing-propeller slipstream or wing-jet interaction. 



3, Conclusions and Recommendations 


The classical formulation of airfoil-jet interaction in inviscid 
incompressible flow has been extended to the linearized inviscid com- 
pressible flow with Mach number nonuni formity. The numerical analysis 
of this new formulation showed the following results. 

(1) The reflection coefficient A.^2 Eq. (12) is the pri- 

mary parameter which determines the aerodynamic characteristics in the 
airfoil-jet interaction in either compressible or incompressible flow. 

A negative value of increases and - ^ • A.n equivalent 

solid body may be used to represent the jet effect only if the flow 
parameters are such that ~ “i* it possible to use 

to correlate the experimental data of wing-jet interaction. 

(2) For the upper-surface blowing configuration assumed in the 
present study, the Mach number nonuniformity with higher jet Mach 
number increased the lift slightly because the reflection coefficient 
becomes slightly more negative compared with the uniform Mach number, 
or the incompressible case. 

(3) High jet temperature has detrimental effect on lift as it 

reduces the jet dynamic pressure, so that positively. 

(4) For the upper-surface blowing configuration, an optimum jet 
thickness exists for maximum lift augmentation at a given thrust. 

In order to determine the effects of airfoil geometry, turbulent 
mixing, non-matching of static pressures, etc., on the results obtained 
in the present study, it is recommended that the finite-difference 





method be used to solve the interaction problem. Because of the 
generality of the finite-difference method, it is expected that the 
solution technique of such general configurations is also applicable 
to determining the Coanda jet performance. The two-dimensional general 
solution, when coupled with a nonlinear lifting line theory, is cap- 
able of providing the three-dimensional aerodynamic characteristics 
for moderate to high aspect-ratio wings. 





APPENDIX C 


A QUASI VORTEX- LATTICE M)2THOD IN THIN WING THEORY 


1. Introduction 

In the conventional vorte:K lattice method (VLM) , the continuous 
vortex distribution representing the wing in a uniform flow is replaced 
by a discrete vortex distribution, where a vortex is placed at each 
elemental quarter chord. Reasonable results have been reported in 
the past [32and35]. However, difficulties exist in the conventional 
VLM: (1) The convergence of solutions (C_ , C and is 

slow with respect to the number of elements used, in particular for 
low aspect-ratio swept wings [32]. (2) The method used to compute 

the induced drag implies that the leading edge thrust is distributed 
over the chord, as has been criticized by Hancock [34j . (3) The pre- 

dicted pressujre distribution is reasonably accurate only away from 
the leading edge. Near the leading edge, the predicted pressure 
level is always too low, even in the two-dimensional case [35] . It 
follows that if the leading edge thrust is computed by using the 
leading-edge pressure distribution, the predicted induced drag would 
be too high. The purpose of the present investigation is to find 
ways through theoretical consideration to revise the conventional 
VLM so that the above-mentioned difficulties can b:^ overcome, 
yet retaining the simplicity of the VLM. 

2, Two-dimensional Theory 

In thin airfoil theory, the downwash on the airfoil is related 
A 

to the vortex distribution through the following equation: 





— L (■' 




U) 


X - x' 


where the coordinates are referred to the chord length, and U/ and / 
are nondimens ionalized with respect to the free stream velocity V^, . 
According to the conventional VLM, Eq. (1) is discretized to be [35] 




W 


where the elements of the matrix [A] is given by 


V 

d - X. 


'/z i- )' - i 


‘ J 

X. = u -- 1) -h 5 ^;t/4 


(i> 




= U' - 0 + '"Vi- 


Assuming that 1> the exact solution of Eq. (1) is then given by 


— z/jHE: 


C^) 


The exact lift and moment coefficients are 


‘^e = T 


— • ^ vi 1 = 2 ( 2 /Ei 

J V... J ^ X 


cix ” 2 7r 


7r 




e.e. 


4- 


C 2 



Let Cj = lim Ui 5 c>y^ . Then the suction force can be computed by 

$ = lcsc^ to 

For a flat plate = 1 . 

To check the accuracy of the VLM, take i = j = 1 in Eq. (3). 

It follows that djj = 2 and hence, /(.\) = TC from Eq, (2) at x = 1/4. 
But from the exact solution, Eq.(4) ^ ^ ^ , so that the VLM 

solution is too low. On the other hand, is exact, because 

It has been shown by James [35] that C_^ and obtained by VLM 

are always exact for any number of elements greater than 2. At least 
two elements are needed to compute . To examine the accuracy 

of the predicted for different number of elements, a few com- 

puter runs were made for 2 through 10 elements. Some comparisons 
are given in Tables la and lb, where N is the number of elements. 


Table la Comparison of Vortex Distribution 

for N = 3 


X 

i by VLM 

Exact / 

0.083333 

5.890486 

6.633250 

0. 416666 

2.356194 

2.366432 

0.750000 

1.178097 

1.154701 




Table lb Comparison of Vortex Distribution 

for N = 10 


X 

i by VLM 

Exact 

0.025 

11.070787 

12.489996 

0.125 

5.244057 

5.291503 

0.225 

3.701687 

3.711843 

0.325 

2.879090 

2.882307 

0.425 

2.325419 

2.326320 

0.525 

1.902616 

1.902380 

0.625 

1.550279 

1.549193 

0.725 

1.233896 

1.231764 

0.825 

0.925422 

0.921132 

0.925 

0.582672 

0.569495 


From Tables la and lb, it is seen that the primary error in the pre- 
dicted Y(p^) occurs near the leading edge and the secondary error 
near the trailing edge. 

Examination of Eqs, (1) and (2) reveals that the Cauchy singu- 
larity and the square-root singularities at the leading and trailing 
edges contained in have been completely ignored in the VLM. 

The conventional VLM was formulated from the standpoint of computing 
the overall aerodynamic characteristics, instead of the local pro- 
perties. Therefore, if the VLM is to be improved, the above-mentioned 

»s 

singularities must be properly accounted for. This being done in the 

A 

following. 

Transforming the x -coordinate to the ^ -coordinate through the 




following relation 


X ^ ( I — cos 0 ^ ^ 


{i) 


Eq. (1) can be written as 




f I ^ ^ 


U?(9) .= - 


Zlt 


(T) 


0 ~ 


Let 


9(e) = (9 

Then, to eliminate the Cauchy singularity in the integrand, Eq. (7) 




is rewritten as 


^ne) = - ^ !3 j e0 . - 3(0) ,,, ^ ( 


3(e) de' 


K Ccs 0 - CCS 9 ' 


-L tKe') - 3(0) 

27T ) 

o CCS G- - CCS 6' 




Now, ^O) does not contain singularities, because the square-root 
singularity in 3"(0) has been eliminated by the factor sin0 . Since 
the integrand in Eq. (9) is finite everywhere, the integral can be 
reduced to a finite sum through the modified trapezoidal rule (mid- 
point rule) [22] . It follows that 
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( ^(e'j — ^(£) jg' ^ Ao' 'Si 
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2 ^ - I 
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f{ = 1 ^ ccii. 9 - Co;. (.ZAjiL 7j-\ 

^ 2t^ J 
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Cos Q — cos 
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where the last step is allowable t£ B ^ 0^ , Now, the primary con- 

cern is to eliminate the last term in Eq. (10) . This can be done 
by choosing a particular set of control points, i.e. ^ - values. 

Let ' zeros of (A) , the Chebychev 

polynomial of the first kind. Then 

U) = K 0\ - A,) ( A - AO ( A - A^,) 

and by logarithmic differentiation, 

^ _ 'T>v'(a) 

r,(A> 

Now, if A. = cos B f it is known that 
T,, ( B 


The zeros of T are therefore 


TT > or 


^ ^ k - I _ 

^ -nr 

A. 


2- 


* = (^TTr = 1 . 2. 


(,;a) 


C , t> 



In order that the sum in Eq. (11) be zero, it is required that 


**T* ^ c( 




^in hid 

— M 


p ' dA 

must vanish. This is true if — i 7T » 


C/.3) 


A • — C-OS 0. -~ 1— - \^= \ 2, - -- -- N I (<>) 

^ I N 

For I = N, T^' ( A^) > -^V^'ces/^Ttand tJ C ) Jt,^ ( A;) ^ - N 

This is the situation when & = 7T in Eq. (6), or A = 1. But at 0 = 7T^ 
3(G) = /(0)i sih 0 — C. Hence, the control points in Eq. (14) can be 
extended to t = N, i.e., the trailing edge. On the other hand, if 
c = 0, i.e., at the leading edge, it can be shown from Eq. (13) that 
/ L,,' — • Since the unKuown function /((?) is known to possess 
a square-root singularity at 0 = 0, it follows that 


\\m C,/IZL .2 — 2d 

d^o 0^0 x-^v ^ ^ 

Hence, if the control points are chosen at N + 1 points, i.e,, = 0, 

1, . . . N, in Eq. (14), (N + 1) unknowns, ^ ^ 0 can be 

computed. Note that the parameter C in Eq. (15) is related to the 
leading-edge suction parameter C^, because » Hence, 

c^ = ii,„ u(.x)/3" = iiw .^ /r = 4;C 

X — > O X — > V ^ 

If the singularity in / is weaker than the square root singularity, 
then (^ = 0, as in the case of parabolic camber. It follows that the 


C.7 



above-described procedures not only predict (jf— values, but also the 
leading-edge suction directly. 

Using the concept just described, Eq. (9) can now be reduced Co 


( 7r i- /ft e. 




^ «=l cciC,. - n'i.0. 
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^^(.cON^ .-.= 0 


4 o 


J_ ^ 

; _ V + 


N(L 
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t = D 


i 0 


Ui) 


where 


t = /(6ft) 


= Tt I - 


-), = ! 2, ■ - ■ - M 

. N' 


Eq. (17) can be written in a matrix form for the solution of ^ as 


follows : 


vrhere 




a. = -!— 

t Ik 

X, = ± C , 0 = 


. N 






. N 


OO 


Cf?) 





Once is obtained, C can be computed as 






The lift coefficient is given by 




— I ^ ^ ^(t>) 5.j>i 

^ s \4.'' ^ 

- 3L ^ / ■ n 


& Up 


where 0^^ is given in Eq. (18). Similarly, the pitching moment co- 
efficient about the leading edge can be shown to be 

M 


c. 


j-— = -zf ,=-l 

i o «<.» . i- 

IL V - 
~TTT 7 “ ^ 


To show the remarkable accuracy of the above formulation, take 
N = 1. Then = 1> ^1^—4: ^^id '/a- • it follows that for 

10“ = 1, = 2. This solution is exact, as can be seen from Eq. 

(4) by putting X = '/^ . The lift coefficient is 




7r 


x-zx j - 2 7r 


and the leading edge suction parameter C is 


Cl - -T i- 4 -a z - a 


(.c) 

(’■i) 

S>V 1 C 



Both are exact ♦ can not be correctly predicted when N = 1. 

For N = 2, the matrix [A] in Eq. (19) can be shown to be 


[A] 


-1 






/ 


The solution for | j j can" be easily obtained to be 

3' = 2 ( It VT) 

^ C/2 -l) 


It follows from Eq. (20) that the suction parameter C is 

d = -L -t -Lj^ 2 ( 2 ( /7 - i/] = 2 


The lift and pitching moment coefficients are 

c^= tiCfT =27t 

^ ( I - r - Crt - 94 ( /t =?)] 

— — 2L 
2 

It can be easily checked that the results for 
are all exact. 

For N > 2, computer runs were made for N 


, C , , and 

up to 10. Again, all 


C. fo 



aerodynamic characteristics of interest are reproduced exactly. 

To check whether the present method will work with the same 
accuracy for cases with camber, consider the case of parabolic camber. 
For 

A(l- X) , 

it can be shown that 




27r 



/a;^( 0 

A'" 



C ^ and C_ are 



Of course, the suction parameter C is zero for this case. Computer 
runs with N = 5 and 10 again reproduced the above results exactly. 

3. Three-dimensional Theory. 

To apply the above theory to a thin wing of finite aspect ratio , 
the expressions for the induced velocity due to vortex segments in 
space are needed. According to p. 43 of Ref. 28, the needed expressions 
are 


(R) - 


^‘7T J 
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r 


(^.5 ) 


for a bounded vortex segment of strength 2^(s)ds of length L and 



for the corresponding trailing vortices. Tlie notations are (refer 
ring to Fig. 1) 


r= I - m: 

y '(0 = 

^ x'-x) ^ ( 3'- y ) i ^ ^ 


= the -coordinate of one end point of the 
bounded vortex element 


\ 




«,-i:) 


' / 
I / 

1 /' 

\;/ 
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Fig. 1 Vortex Segment Geometry 



Eq. (23) can be integrated over the bounded vortex segment L to give 


(See Appendix A) 




ib'i 



where 


+-^(^1-^^*’^ 

t'= ix,-x)t -t- (%.-%) i 




N) 


Eor the application in thin wing theory, it can be shown that 

a'xl' := pi ^ - u.-x,)( 

i>' • t = - x) 

It follows that Eq. (25) becomes 










( X,) ( X, - X) -f /x\ ^,- ■'j , X 


(^7) 


Assuming straight trailing vortices on the wing plane, the integral 


in Eq. (24) can be reduced to 






r i- - 1 




) 



Hence, the downwash due to the trailing vortex from the first end 
point ( ) of the bounded vortex segment is 




^•?T 















In the same manner, the downwash due to the trailing vortex from the 
second end point ( of the bounded vortex segment is (noting 

that f. = c ) 




f 

Tir 


), a,-!) 

Xm 






Thus, the total do\mwash produced by a skewed horseshoe vortex element 
is the sum of U,'; , • 

For a rectangular horseshoe vortex element, and Eq. (27) 

shows the usual Cauchy singularity in the x direction. For skewed 
elements, with ^ inside the element under consideration, and if ^ 

(the y-coordinate of a control point) is chosen so that 

y,- = k'(^,- so 

where K**is some constant less than one, the expression in the denomi- 
nator of Eq. (27) can be written 

= X, tKXx..- X,) -x] 

= '>0 


C, '1- 



which vanishes on the vortex segment at the chosen y— coordinate . Thus, 
the Cauchy singularity exists and the method of integration described 
in the two-dimensional theory is applicable. 

For notational convenience, let 

/ (?., - >-)> - 'if 

u, U'. 3) = r 1 - — 1 

j) ,= I I - .. „ 1 (3<^) 

^=■-3- 


The total do^-mwash due to a strip of vortex distribution is therefore 
given by 


Ur(x J- ( 

^ '^'ir ) 
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(36-) 


Transforming the integration variable to a non-dimensional one based 
on the local chord and Chen to the 0 -coordinate, Eq. (35) can be 
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reduced to a finite sum as follows: 






4-tt a/ 
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where 


X. = 










\ = 


± 


T ) CAJ 3 ( 

1 

i 

L * \ 



)] , ^ = 


^ — the leading edge x~coordinat at y^ 


V _ 

the leading edge x-coordinate at 




= chord length at y, 


— chord length at y^ 


So far, the control-point locations have not been chosen yet. In 
accordance with the consideration in the two-dimensional theory, the 
x-coordinates of the control points should be chosen so that 


C. t 6 







f-Li- 


co,s 




C —\ 


f\l 


UJ ) 


where is the leading edge x-coordinate on the chord C through 

the control point. To choose the y control location, it should be 
noted that both Eqs. (33) and (34) involve the Cauchy-type singular- 
ity in the y— coordinate. Since the downwash expression given in 
Eq. (35) is to be summed over the spanwise strips and the summation 
is equivalent to integration in some sense, the spanwise summation 
can be more accurately performed if the singularity can be removed. 
For example, the summation 


■um ' 




can be accurately performed, where hi ^ ) is either of the 


expressions in the brackets in Eqs . (33) and (34). Now, if ^ ^ ^ 


1 > 






i ^ - 'i 


- ~ t, a) ^ 


j ij-'i 


and the last term can be eliminated if 


J a. a, y 
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The choice of the control points and the spanwlse strips is illus- 
trated in Fig. 2 for M = 5, Admittedly, this choice does not elimin- 
ate the effect of Cauchy singularity, because both trailing legs y^^ 
and y^ do not occupy the left side of the first strip and the right 
side of the last strip simultaneously. However, it should reduce 
the error, as has been proved in the numerical study. 



Fig. 2 Illustration of Control Points and Spanwise Strips 

To find the leading edge suction parameter, note that the leading 
edge suction term comes only from the bounded vortex distribution, and 
hence, from Eq. (27). Transforming Eq. (27) to the Or -coordinate, 
and noting that Eq. (31) can be written 

(X,- 3 ,)- 3 ) = ( 3 ,- 3 , (005.6 


C. IS 



since 


^. + k'(x^-a:,) x,-i- \t-7( I - 

X — 4 ^ ( I _ 0 ) 

c. > 

it is found that 




U,(e^>/(0*; ^ s\y,d' ^ 0 ' 

( vcs9 


- -L j£^ S S 

^ ^=il ( v^.) (c^s e - 


{M) 


in accordance with the concept discussed previously in relation to 
Eqs. (9) and (10) ♦ 

As shown before. 
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it follows that 
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Thus, the second term in Eq. (4T) becomes 


I im !— ^ ^t( 0) ifiV) ^<n & 

^•*>0 ^_, ( H.) (^C-6.S 0 -CAiS0^) 
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If Eq, (45) is used in Eq. (36) which must be summed spanwise, it 
is found that 


N/<i/ t A* 




+ ^ ^ 

^ j- S'" 

^ ^ - 0< ) . a+ A = ;i<, ^"^9 
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where the unmarked summation is to be performed in the spanwise dir- 
ection. Once tlie parameter C for each strip has been obtained, the 
magnitude of the sectional leading-edge thrust coefficient can be 

computed [36] : 


C. 


2l 

8c 




y j - 




where 


l.-», AC„(x-X^f = 1:„ z/u; =2CIZ 


Hence , 

K 

For a skeined bounded element, the total vortex, density of the 
element is /^^/cos'Y, where 4^^ is the strip width and is the 
element sweep angle. Its lift-producing component is obtained by 
multiplying by cos . It follows that the lift-producing vortex 
density per unit span is ^ . The sectional lift coefficient is 
therefore 

= T'Ta . \ 

^ If ^ y 

~ ^ J ^ ^0 - —j— ^ s.y. 


C.^l 



Similarly, the sectional pitching moment coefficient referred to the 
y-axis is 




where ^ was given in Eq. (37). The sectional drag coefficient is 
= Cj- ■ o( — 

To compute the total thrust coefficient, it is assumed that the 
predicted represents the local value, instead of the average 
strip property. The total value will be computed by integration: 



where the integration Kaa been reduced to a finite sura by the conven- 
tional trapezoidal rule and 



c - I , 


M - 1 




Note that in the conventional trapezoidal rule, the values of the 
integrand at both integration limits appear in the finite sum. In 
the present case, both values are zero, because sin = 0 at <|> = 0 
and . The values of 4^. in Eq. (52) coincide with the y-control 

locations given in Eq. (40) . It follows that the total thrust coef- 
ficient is given by 
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similarly, the total lift and pitching moment coefficients are 
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To check the accuracy of the above formulation, four planforms 
were chosen. The convergence of solutions for the rectangular plan- 
form of aspect ratio 2 is illustrated in Fig* 3. It can be seen 
that results with 3~digit accuracy can be easily obtained with the 
present method. Comparison of aerodynamic characteristics predicted 
by some analytical methods is made in Table 2. The values of 
and given by the - • > 


Table 2 Comparison of Aerodynamic Characteristics 


of Rectangular Wing of A = 2 by Different Methods at = 0 



rad'* 

rad"' 

.. . 1 

Sii/Cl (Near Field) | 

_ . 1 

Present 

N = 6, N = 15 
c s 

2.4751 

-0.5185 

0.1591 

Present 

N = 5, N = 20 
c s 

2.4744 

-0.5185 

0.1589 

VLM 

N = 6, N = 20 
c s 

(Ref ,32) 

2.5239 

-0.5334 

0.1554 

NLR Method 
(Ref . 3?) 

2.4744 

-0.5182 

■ i 

0.1609 


Note; N = chordwise number of vortices. N = number of spanwise 
c s 

vortex strips^ = 0.15915 . 





NLR method are regarded to be reliable, because the methods of BAG 
and NPL produced the same values [37] . It is seen from Table 2 that 
the present method gives more accurate results than the conventional 
VLM. 

The second planform is the Warren 12 planform which was also 
evaluated in Ref .f37] (see Table 10 of this reference) . Table 3 shows 
the results. Again, the present method seems to give results closer 


Table 3 Comparison of Aerodynamic Characteristics 
of Warren 12 planform of A = 2^1 by Different 


Methods at Mc©= 0 



C^ rad"' 

Lot » 

C„ ,rad"’ 

^Di^C^ (Near Field) 
L 

Present 

Nc = 5, Ns =10 

2.7515 

-3.0936 


Present 

Wc “ 5, =15 

2.7468 

-3.0898 


VLM 

Nc = 5, Nc, =15 
(Ref .32) 

2.7944 

-3.1775 

0.1055 

NLR Method 
(m = 15) 
(Ref .37) 

2.7373 

-3.1074 

0.1201 

NLR Method 
(m = 31) 
(Ref .37) 

2.7576 

-3.1155 

0.1135 

NPL Method 
(m = 15) 

(Ref .37) 

i 

2.7270 

- 

-3.1038 

0.1227 


Note: m = Humber of collocation sections over the whole wing. 

l/(TrA)= 0.1125 


C. 25 





























to those given by the NLR method. 

Another low aspect-ratio planform evaluated is the 45®-sweep 
constant-chord wing of A = 2. This planform has also been investi- 
gated in Ref. 32. The characteristics of solution convergence are 
illustrated in Table 4. 

Table 4 Aerod 3 maTnic Characteristics of 45®-sweep 


Constant— Chord Wing of A = 2 at = 0 , 4 


Present 

Method 




VLM (Ref. 32) 

Ns 

10 

15 

20 

25 

25 

C- rad 

MX » 

2.2658 

2.2604 

2.2583 

2.2573 

2.2924 

C^ rad ■’ 

-1.5136 

1 

-1.5106 

-1.5096 

-1.5090 

-1.5468 

Wc* 

(Near Field) 

0.1525 

0.1556 

1 

0,1573 

0.1583 

0.1520 


Note: '/(TtA) = 0.15915 


Again, 3-digit accuracy for C and C^ is easily attained. However, 

Lo< ^ 

the convergence of the near-field z seems to be slower than 

L 

the other planforms evaluated. 

To assess the accuracy of the method for high aspect-ratio wings, 
the rectangular wing of A = 7 is chosen. The aerodynamic characteris- 
tics of such a wing are shown in Table 5. It is seen that the results 
change little when the number of strips is increased from 10 to 15 


with N ^ = 5 . 
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Table 5 Aerodynamic Characteristics of 
Rectangular Wing of A = 7 at = 0 . N^. = 5 



Present Method 

VLM (Ref. 3 2) 

Ns 

10 

15 

15 

C rad 

* 

4.4235 

4.4215 

4.5103 

rad 

® > 

-1.0647 

-1.0642 

-1.0888 

Si/C‘ 

(Near Field) 

0.04647 

0.04646 

0.045 


Note : 0.0455 


From the above evaluation, it may be concluded that the present 
method seems to yield consistently the results accurate to 3 digits 
for C and , within the range of number of strips used in the 

investigation. For low aspect-ratio wings, in particular with sweep, 
the convergence characteristics for the near-field with re- 

spect to the number of elements have been improved over the conven- 
tional VLM. However, higher accuracy seems not easy to attain with 
reasonable number of elements. It is thought that this may be due 
to the highly three-dimensional nature of the flow field associated 
with swept low aspect-ratio wings, so that the assumption of con- 
stant vortex strength over a strip is less applicable. 


CM 




In order to treat complex planforms, such as variable-sweep 
wings, and wings with rapid variation in section properties, such as 
wings immersed in propeller slipstreams, it is convenient to intro- 
duce the following scheme of spanwise vortex strip distribution. Any 
spanwise interval on a semi-span may be divided into vortex strips 
by the semi-circle method described above. This has the advantage 
that the region of rapid variation of section properties can be covered 
with finer vortex strips. This is illustrated in Figure 4. Gaps 
between spanwise sections resulted from this procedure must be eliminated. 



Figure 4. Second Scheme of Spanwise Vortex Strip Distribution 

Any interval can be transformed into [o,7r] on the 4?-plane, For 
example, any station y in the interval expressed in 

terms of ^ through the relation: 

CO.*, 

' Z 2 ' 

Thus is transformed into [Oj'n']. Let the vortex strips be 



chosen so that 


y, - -g. 

fc/z -i, 




^3 - I 




M, 


C^-7) 


and the control points given by 


Vc ~ I 

Wa-“ 



cos 


cTr 



‘ = I, 


M,- 




SliDllar consideration is also applicable to the interval [o,^l] with 
M 2 replaced by It follows that the total lift coefficient can 

be computed as follows; 





c U y 
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Si^i<j:> d<j-> -j" 


z 
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The total pitching moment coefficient and the total thrust coefficient 
Cx can be found in a similar manner. 

In the thin wing theory, the convergence of solutions is indicated 
by the agreement of the near-field and the far-dield induced drags. 

The far-field induced drag is given by [Eq.l43 of (38)]: 



where is the induced angle of attack given by [Eq. 61 of (.30)] 


•- STr J, 


c,.c 


iXy* 


-y. 


C/-y'7 


Let 


y'-yV . j'=-c-oi.f 


(£i) 


Eq.(61) becomes 


o(. = - ^ 

<i-vr 


iCi) 


■I '■'i-vy 

To Integrate Eq.(63), Multhopp's interpolation formula is used for c^c [38] 
Thus, 


M-1 






C^<^) 




where M-1 is the number of interpolation stations over the whole span, 
and 

<h — — 2E , » — 

/V| ' , / • I 

Substituting Eq.(64) into Eq. (63) , it is obtained that 
cp(. =r -- — - — ^ (c, c) s;«a4; \ 

“i-Ti-b M " A-I 3—: y;i 

' O ^ — c-<, s 4> J 


(^C5) 


:= _! V ^ 




h=j 


S f'-n <i> A=| 


>v. sin a 


iU) 


where the following relation 

c'tr 


I siv\a_4* 4*^ a4* 


(f^T) 


has been used (Appendix IV of Ref. (38)). Now, choose <|> such that 


C.^0 



<b 

" M 


It follows that (Appendix IV of Ref, (30)) 

^ (^l)' 

3ir^ A<t = — 


' I 2 Sin c}:^^ 

M ( c<5i. ^ ._ coii ^ 


Eq»(66) becomes finally 


/ V i . 


OncecXi's are obtained from Eq,(70), Eq,(60) can be integrated 
for the induced drag, following the scheme of Eq. (53) . Thus, 

^ ^ ^ ( d.e':- o<c)^ 

where is given by Eq. (68) with V replaced by i. 

In Eq, (70) , 's and c’s at stations specified by Eq. (65) are 
required. However, only those at stations given by Eq. (58) are known. 
Therefore, some interpolation procedure must be used. In the present 
method, the following Lagrange interpolation formula was used: 

C - Cy-y.)(y-/^)---(y-y„) ^ (y-y„) ^ ^ 


with similar formula for the local chord c . 


In the following illustration of some computed results, M is 
taken to be 31, except for the variable-sweep wing where M = Al. It 



Is assumed that the convergence of solutions is achieved if the near- 
field induced drag differs from the far-field value by less than 1%. 

The rectangular iving of A = 2 is again employed to make the comparison. 
This is illustrated in Table 6* The second planform chosen for com- 
parison is that of a delta wing of A = 2 at 0.13. The comparison 

is made in Table 7 with experimental results and other analytical 
methods. It should be noted that the vortex arrangement for the VLM 
is judged, from the study made in Ref.(3E), to be the best for a total 
of 120 vortices. The near-field induced drag parameter by the VLM is 
3.45% too low. The present results are comparable in accuracy with 
Wagner’s method. However, Wagner’s program requires large ammount of 
processing time, being about 8 times higher than the present method 
( 71.6 sec. in the Honeywell 635) for this configuration. 

The third planform is the sweptback and tapered planform of A = 8.02 
given in Ref, (39), The geometry is given in Fig. 5. The results are 



Figure 5. Geometry of Planforms for Tables 8 and 10. 
shown in Table 8. Additional results for other planforms are compared 
in Tables 9, 10 and 11. In all cases, reasonably good results are 


C . 



obtained by the present method. One important feature of the present 
computer program is that Purchell’s method [21] has been used for 
solving the matrix equation, instead of the standard matrix inversion 
method. This requires only 17K (decimal) memory and 95.8 sec. for 
execution in the Honeywell 635 for a total of 120 vortices, compared 
with 26K (decimal) memory with overlays and 350 sec, processing time 
for the VLM of Ref, 32. A fraction of the latter processing time may 
be due to the computation of additional quantities, such as the side- 
wash, etc. 

The spanwise distribution of the sectional lift coefficient and 
the sectional induced drag coefficient is illustrated in Fig. 6 for the 


variable-sweep wing. 



Table 6* Aerodynamic Characteristics of a 
Rectangular Wing of A = 2 at = 0 


Method 

Lo( ,rad 

C , rad ^ 
moc 

S../Cl2 
11 ^ 

\i/% 

Present 

Nc=8,Ns=15 

2.4707 

-0.5173 

0.1595 

0.9972 

VLM 

N^=6,N„=20 
(Ref. 32) 

2.5239 

^ — 1 

-0.5334 

0.1554 

0.9747 

NLR 

(Ref. 37) 

2.4744 

-0.5182 

0.1609 

1.0101 

Wagner 
(Ref .40) 

2.4778 

-0.5180 

0.1619 

1.0167 


Table 7. Aerodynamic Characteristics of a 
Delta Wing of A = 2 at = 0»13 and <x = 4.3' 


Method 


Present 

Nc=3,Ng=35 


Lamar 
(Ref .39) 


VLM 

N^=6,N.=20 

0 . 1 ) 

(Ref .32) 



0.1649 


0.168 


0.1654 


■k 

a. c . 

S. ./C 2 

ii' L 


0.3767 

0.1625 

1.0078 


0.363 


0.3823 



0.1549 


0.9655 


Wagner 
(Ref .40) 


0.1663 


0.3783 


0.1576 


0.9795 


Experiment ^ 

(Ref. 59) 


0.396 



measured from the leading edge of the mean geometric chord as defined 
in Ref. 39 and referenced to the mean geometric chord. 



















































Table 8. Aerodynamic Characteristics of a Sweptback and 
Tapered Wing of A = 8.02 at = 0.19 and o< = 4.7® 



See the footnote in Table 7. 


Table 9. Aerodynamic Characteristics of a 
Rectangular Wing of A - 7 at 0.0 


Method 

, rad 

C .-1 

m^,rad 

11 L 

Si, /Si 

Present 

N =5,N -15 
c ’ s 

4.4142 

-1.0624 

0.04658 

0.9972 

VLM 

N^=5 ,Nc;=15 
(Ref. 32) 

4.5103 

-1.0888 

0.0450 

0.9597 

Wagner 
(Ref. 40) 

4.4300 

-1.0680 

0.04888 

1.05159 


Lamar 
(Ref. 39) 


4.40 


-1.047 
















































Table 10. Aerodynatnic Characteristics of 
Warren 12 Planform of A = 2^2 at Mro = 0 




Method 

, rad 

C .-1 

m^,rad 

/C ^ 


Present 

N^=5,Ng=20 

2.7382 

-3.0844 

0.11428 

1.0054 

VLM 

N^=5,Ns=15 
(Ref. 32) 

2.7944 

-3.1775 

0.1055 

0.9275 

NLR 

(m = 15) 
(Ref .37) 

2.7373 

-3.1074 

0.1201 

1.0564 

NLR 

(m = 31) 
(Ref .37) 

2.7576 



-3.1155 

0.1135 

0.9921 


Table 11. Aerodynamic Characteristics of a 
Variable-Sweep Wing of A = 4.303 at Mco~ 0*23 



See the footnote in Table 7 




















































APPENDIX D 


Expressions for and Gj^ 


1. Wing-Slipstream Interaction 

According to Eq.(4.33), is given by 



; r 

I Pi 



( 1 ) 


The circumference of the circular slipstream is divided into a 
number of equal straight line segments, so that the circular slipstream 
is approximated by an inscribed polygon. The control points (x, y, 2 ) 
are chosen at the midpoints of the elements. In terms of the cylindrical 
polar coordinate system, as sho^vn in Figure 1, the endpoints of any 
element, (x-j^, y^^, z^) and X2» expressed as follows: 

XI = X5 

Xl = + Pcos (j)j^ (2a) 

Zf + r sin (j)j^ 

^2 = 

V 2 = + rcos (j >2 (2b) 

Z 2 = + r sin ^2 

where (x^ , y^ , z^) is the center of the circular slipstream at the section 
under consideration. With the geometry of the vortex elements defined, 
the various terms in Eq. (1) can be computed in the rectangular coordinate 
system as follows: 



(3a) 


axjt -i 1 yX2^-a,) 

t I [(*,-. aX 9,) 

|ft.'x^'| =L fi" { [ (*»- ^.)(l»rs) - (A|-a)( v-’J,)3 + [(^.'^i)(*^'~^^ 
-(*,-»)( ^’■|:( »^' OiX*-.- *)- (i«- 

b'.? =• C)i^-xX^*-;c,)-i-yS'"('j.-^)(3.-3,).+ 
t'.r = (a:, ,- bX X'S,) t /s"(i,-2-)(».-*.) 
itl = >X+ /3‘ ( )).- >/+ /SX - ^x} (3e) 

la'I = t (a, ( 3f) 



Figure 1. Jet Polar CoorcJinate System 


According the Figure 1, the unit normal vector is 

A ^ '*•. 

=. 5 «|> i- S|n (4 ) 

It follows that 

-)A -l' I ‘^^ 4 ' 


a,)3 


iP| 


■-iMi 


V 

I 


a'KV 


(5) 


In the present linear theory, the slipstream axis is assumed to lie 

^ A 

along the x-axis. Therefore, = i. It follows that 







S I 

i I P| li^'i J 


I A'XX'I 

To evaluate Gj^, it is assumed tliat the trailing vortices are 
straight lines along the x-dirction. It follows from Eq.(4.34) that, 
for trailing vortices emanating from the second endpoints. 


( 6 ) 




^3 






f ^ - c n'- 


i j R’ 




(7) 


where it has been assumed that each trailing vortex line is parallel to 
the x-axis at a constant y. The integral in the last expression can be 
integrated as follows; 


S ^ -5 


d?c' 




7 . 5 - 








>3 



Hence, 




e, - a 




5 


- a 






:1 


X' - a 


(9) 




For trailing vortices from the first endpoints, it should be noted 
that d^i = -d‘C 2 . It follows that Eq. (9) is still applicable in this 
case except that the subscript 2 is to be replaced by 1 and there is a 
sign change. 

2 V . X, -X 






j ]-- ^ 




.i 


+ ft ■-, \ I - : 


X, ~x 


Adding Eq. (9) and (10) gives the final results for G^: 


( 10 ) 




U-^) 


( 11 ) 


• n and G^ • e^ can be obtained by simply taking dot products of (H) Witi 


A» -*v 


Eq.(4) and i, respectively. Since e^ = i, obviously Gj^j^ • e^ = 0. 

A ^ 

To evaluate the downwash on the wing surface, note that n^ = k. 






— - 




^ lb' I lA'l ^ 


(12) 


~ ? y I ^ \ 


Vil 


f -j ii_2 y , _ 3 


( 13 ) 


X, -a 


Hence, 



2. Upper Surface Blowing Jet-Wing Interaction 

Eqs.(12) and (13) are still applicable in this application, 
because the jet has been assumed to be flat with constant thickness. 
However, for a rectangular jet with the control points on the sidewalls, 
n would be j on the outboard sidewall and -j on the inboard sidewall. 

It follows that on the sidewalls 


M • n C .. ) (1^) 

I i'l ^ IbM I ^'1 ^ 

• -n = + — — ^ J \ — — 


(15) 


- 2 


tJ'- 


X, -X 




where the upper sign is for the outboard sidewall, • e^ is still given 

by Eq,(6) and =0, 



APPENDIX E 


Numerical Accuracy for Evaluating 
the u-Induced Velocity 


According to Eq.(6.28), the u-induced velocity due to a particular 
strip of vortex distribution can be accurately evaluated by rewriting the 
expression as 




( 1 ) 


There is no difficulty in evaluating the first integral numerically 
because the integrand is at most Cauchy singular. In the present computer 
program, the second integral is approximated by a finite sum as follows 
by setting x' -b c ( i — <^'^90/2 ; 


a'’ i^_ia 

F(x, j I ax' = Y ^ 


(2) 


To check the accuracy of Eq.(2), consider a rectangular wing with unit 
chord length. In this case, the bound vortex is unswept so that 
Xj^ = X 2 = x’ . According to Appendix D, it can be shown that 


M — Nf ^ “ C = 


.'e = S 






t (X, - 3r 

|«'xrr=^' 1 




= 'i.- 5, 


El 



It follows that 


_ . r (■' 


= 

' ^7r 


^ ( -. \ 'i- 

O <- * r f ,;V« 






I 4x' 

= ± 5 W ' ^ ^1 

h t/» ( av- s/j ‘ * [ (x'-xr/tys‘( 3^- ‘/‘ 


o 

= { W c^^(a.-a) C3.-3) 




■Ii.r^ )^‘^r ifl 


T 


A_Oii_ “■?) 


^ C ( I c ^ - g ^ t N-/s^c ^ ^ f -j. ^ \ 

(3) 


where the fact that 72 - y > 0 and yj^ ^ y < 0 has been used. To check 
the accuracy of the approximation by Eq. (2) against Eq.(3), let yj^ = 1, 

72 = 1.^, y = 1.2 and 3=1, For 5 control points (i.e., M = 5) and 
40 integration points (i.e., N = 40), the results are generated and 
compared in Table 1. It is seen that the present interdigitated integra- 
tion scheme is sufficiently accurate. 



Table 1. 

Accuracy 

of Interdigitated 

Integration Scheme 


yi = 1, 

Y2 = y = 1*2, 

3=1 

X 

z 

F (x,z) by Eq. (2) 

Exact F(x,z) by Eq.(3) 

0.09549 


0.21551 

0.21543 

0.34549 

0.15 

0.27348 

0.27347 

0.65451 

0.27348 

0.27347 

0.90451 


0.21551 

0.21543 

1.0 


0.14537 

0.14526 

0.09549 


0.37473 

0.37465 

0.34549 


0,41389 

0.41411 

0.65451 

0.05 

0.41389 

0.41411 

0.90451 


0.37473 

0.37465 

1.0 


0.21062 

0.21022 

0.09549 


0.43349 

0.43452 

0.34549 

0.025 

0.44171 

0.45643 

0.65451 

0.44171 

0.45643 

0.90451 


0.43349 

0.43452 

1.0 


0.23063 

0.22981 
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